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XXI. Magnetic Hysteresis in Superconducting Colloids 


By A. B. Preparp 
Royal Society Mond Laboratory, Cambridge* 


[Received November 2, 1951] 


ABSTRACT 


The irreversible effects exhibited in the magnetization cycle of a 
superconducting colloid are analysed in terms of a simple model of a 
superconductor. Expressions are derived relating various critical field 
strengths to temperature and radius of the particle ; there is fair agreement 
between experiment and theory. 


§1 
SUPERCONDUCTORS display under different conditions a somewhat 
bewildering variety of irreversible effects when subjected to cycles ‘of 
magnetization and demagnetization. These may be classified into two 
main groups, those in which the Meissner effect is exhibited (i.e. complete 
expulsion of magnetic flux from the superconducting phase) and those in 
which it is not, and it is with the former group, the types of irreversibility 
characteristic of pure superconductors, that we shall be concerned here. 
These may be most simply described with reference to a long rod, of 
effectively zero demagnetizing coefficient, in a magnetic field parallel to 
its axis. Under ideal conditions, if the magnetic field H is increased from 
zero the rod remains superconducting, with the induction B=0, until the 
critical field H, is attained, when it becomes normal and thereafter B=H ; 
on lowering the field the Meissner effect occurs at the critical field and B 
once more becomes zero. The process of magnetization is then completely 
reversible, but in practice this ideal behaviour is rarely observed. The 
most striking, and most easily demonstrated, difference between real and 
ideal behaviour is the effect usually referred to as ‘ supercooling ’ ; if the 
field is lowered slowly from a value greater than H,, the normal state persists 
to fields less then H,, and it is only at some lower field strength, which is 
critically dependent on the nature of the specimen, that the transition to 
the superconducting state occurs. The converse phenomenon of © super- 
heating ’, i.e. the persistence of the superconducting state in fields greater 
than H,, is more difficult to demonstrate, and it is only recently that it 
has been shown unambiguously to occur in macroscopic samples (Garfunkel 
1951). Ifthe sample is not in the form of a long cylinder and has a non- 
vanishing demagnetizing coefficient the ideal reversible transition begins at 
a field less than H, and proceeds through the intermediate state over a 
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range of field strengths. It is now more difficult to discern whether 
superheating occurs, but it is probably observable in fairly thin cylinders 
in a transverse field under favourable conditions (Pippard 1950). On the 
other hand supercooling is still readily detectable, and is almost invariably 
observed in experiments on pure superconductors (e.g. Désirant and 
Shoenberg 1948). 

If we turn now to much smaller samples, such as spherical colloidal 
particles, we find under certain circumstances much more pronounced 
irreversible effects, as is shown by Shoenberg’s (1940) study of mercury 
colloids. The larger particles, of radius about 10-4 em, show considerable 
supercooling, and probably superheating as well, for the area under the 
magnetization curve is somewhat larger than is to be expected (Pippard 
1951). On the other hand, no hysteresis was observed with much smaller 
particles, of radius probably little more than 10~® cm. 

There is not really enough experimental evidence to allow the size- 
dependence of the irreversible effects to be described with certainty, but so 
far as can be seen the general features appear to be as follows. Spheres of 
radius less than about 5 x 10~® em show very little hysteresis ; as the size is 
increased hysteresis appears and becomes more and more marked until the 
radius is a few times 10~-° cm ; in this range of sizes the hysteresis is less 
near the transition temperature than at lower temperatures ; for larger 
particles and macroscopic samples the hysteresis is less, and at the same 
time becomes smaller at low temperatures than near the transition temper- 
ature. In the following sections. we shall investigate the extent to which 
these observations are consistent with current views of the nature of the 
superconducting state. 

§2 

The extensive study by Faber (1952) of supercooling in macroscopic 
superconductors has elucidated the general features of the way in which the 
supercooled normal state breaks down. The superconducting phase is 
always initiated at one or more centres, where apparently superconducting 
nuclei are more readily formed, and the existence of super cooling argues 
the existence of a potential barrier which must be overcome before the 
nucleus can be created. Once a nucleus has been formed it may expand 
without further hindrance to convert the whole sample into the super- 
conducting state. According to Faber the surface energy at the interface 
of normal and superconducting regions is largely responsible for the 
potential barrier, and his observation that super cooling is most pronounced 
near the transition temperature is consistent with the temperature 
variation of the surface energy, as deduced from studies of the intermediate 
state in small superconductors (Désirant and Shoenberg 1948). 

Now it is likely (Ginsburg and Landau 1950; Pippard 1951) that this 
interphase surface energy is a consequence of the impossibility of a sharp 
phase-boundary ; in pure metals the transition layer between the two 
phases may be as wide as 10-4em. According to the views of the present 
author (Pippard 1951) this is a consequence of the long-range order 
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exhibited by superconductors, which may be expressed in the following 
way. If we suppose the superconducting phase to be characterized by an 
internal parameter w, the degree of order, which takes a non-vanishing 
value in the superconducting phase and is zero in the normal phase, then 
# cannot vary at will from point to point, but only slowly, so that, for 
example, at a phase boundary in a pure superconductor it cannot drop to 
zero in a distance less than about 10-4 em. 

On this view it is clear that in a colloidal particle of diameter less than 
10-4 cm the transition from one phase to the other cannot occur by the 
creation of a phase boundary, but must affect the whole particle at once. 
There is no very definite experimental evidence on this point since the 
measurements have necessarily been performed with an enormous number 
of particles simultaneously, and on account of the size-variation of the 
particles, if for no other reason, the measured transitions into the normal 
state are spread out. However the series of measurements by Shoenberg 
(1940) on mercury colloids of progressively greater uniformity lend some 
support to the hypothesis that the transition of a single particle is accomp- 
lished suddenly, without the appearance of the intermediate state. We 
therefore interpret the change in hysteresis behaviour which sets in at a 
particle size around 10-4 cm as associated with the possibility in larger 
particles of an appreciable variation of w over the particle, or even the 
creation of a phase-boundary. No attempt will be made here to discuss 
magnetic hysteresis in the larger particles, and we shall confine the dis- 
cussion to particles sufficiently small that it is reasonable to suppose w 
to take the same value at all points. 


§3 

It is not difficult to see how hysteresis arises in the magnetic behaviour 
of these colloidal particles. In the absence of a magnetic field the Gibbs 
function, G, (which in this case is the same as the Helmholtz free energy) 
of such a particle may be supposed to depend, among other parameters, 
on w and the temperature. For example, the dependence on w at some 
given temperature may take the form of curve (qa) in fig. 1, and the equili- 
brium value of w is that corresponding to the minimum of G. Now ina 
magnetic field there will be a positive magnetic contribution to G, which 
will depend on the extent to which the magnetic field is excluded from the 
particle. Since we may expect that the penetration depth, A, will be 
closely linked with w, becoming infinite as @ tends to zero, the magnetic 
contribution to @ will vary with w somewhat as shown in curve (6) ; the 
larger the particle the steeper will be the initial rise, and the curve will 
flatten off as A becomes small relative to the size of the particle. The 
- Gibbs function itself is the sum of the two contributions, as shown in curve 
(c). The curves have been drawn to represent the state of affairs for such 
an applied field that the minimum of curve (c) at w’ hae oe same value as 
at the origin, that is, the superconducting state (w’A0) is in thermo- 
dynamical equilibrium with the normal state (w=0). However, under 


U2 
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these conditions no transition between the states can occur, if we ignore the 
possibility of thermodynamical fluctuations, on account of the maximum 
of G between the minimum and the origin. In order to allow a transition 
from the superconducting into the normal state, the magnetic field must be 
increased, i.e. curve (b) must be scaled up in proportion to H?, until the 
potential barrier disappears. Conversely, in order to allow a transition 
from the normal into the superconducting state the magnetic field must be 
reduced below the value corresponding to thermodynamical equilibrium. 
Thus under the conditions for which the diagrams have been drawn we 
should expect both superheating and supercooling to occur. This is 
exhibited diagrammatically in fig. 2, which shows the variation of G with 


Fig. 1 


(6) 


O A) | 


Gibbs function of small particle : 
(a) in absence of field, 
(6) field contribution, 
(c) resultant curve for G. 


w for various values of H. The top curve is that which just allows a 
transition from superconducting to normal, and the bottom curve that 
which just allows the transition from normal to superconducting. From 
this diagram we can construct a curve showing how w’, the equilibrium 
value of w in the superconducting state, varies with H. Such a curve is 
shown in fig. 3 in which the arrows distinguish the branches corresponding 
to increasing and decreasing fields. It is clear from this figure that there 
is a discontinuity in w at both the superheated and the supercooled | 
transitions. ; 

If the particle is rather smaller, having dimensions comparable to 
Xo, the penetration depth at 0°x, the field is only partially excluded for all 
values of w, and the magnetic contribution to G follows a curve which 
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sae much less steeply than fig. 1 (b). Under these conditions there may 
be no potential barrier to give rise to superheating and supercooling. This is 
illustrated in figs. 4 and 5, which are analogous to figs. 2 and 3, respectively. 


Fig. 2 


O ; &) | 


G versus w for different values of H. The numbers attached to each curve are 
relative values of H. 


Fig. 3 
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H (arbitrary units) 
Variation of w’ with H, from fig. 2. 


The equilibrium value w’ now is a continuous function of H, and no 
hysteresis is to be expected. Thus at any given temperature we should 
expect to find a critical size of particle below which no hysteresis is 


observed. 
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§4 
After this qualitative discussion it is a simple matter to formulate the 
problem quantitatively. We shall suppose that for a given colloidal 
particle at a given temperature the free energy in zero field is a function of 
w only ; ie. F=F(w), and that the magnetic susceptibility of the particle 


Fig. 4 


e) (23) | 


G versus w for different values of H. The particle size is smaller than in 
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_ _ (arbitrary units) 
Variation of w’ with H, from fig. 4. 


is also dependent only on w. We may therefore write for the¥magnetic 
moment in a field H, M=—2K(w)H. The Gibbs function, G(=F—MA) 
then takes the form 


3 


G(w,H)=F(w)+ HK (w) . Spe Seay 
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As stated before, w’ is determined by the condition (dG/dw),—0. We 
shall now derive expressions, in terms of F and K , for the three fields of 
particular interest, H, the field at which the superconducting and normal 
phases are in equilibrium with one another from a thermodynamical 
point of view, H,, the field at which superheating ceases, and H ge the field 
at which supercooling ceases. 

The field H, is determined by the condition 


G(w’, H,)=G(0, H,)=F(0), since K(0)=0. 


Now it is clear from fig. 2 that if H is less than H, it is possible to find a 
value of w, not the equilibrium value w’, such that the equation 


(Gra EL) =30( 0), ok, ei Wet at cotta) 


is satisfied. On the other hand, if H>AH,, no solution of the equation is 
possible, since the minimum at w’ corresponds to a higher value of G than 
F(0). Thus H, is the largest value of H for which a solution of (2) is 
possible, and consequently, from (1), H,2 is the maximum value that can 
be taken by the function [F'(0)—F(w)]/K(w). We may write this result in 
the form 


H2= Max { (3) 


0<o<1 


BUFO} 
K(w) 


If the analytical forms of / and K are known it is a simple matter to 
compute H, from this formula. 

A similar expression for H,, may be obtained immediately, since it is 
clear from fig. 2 that Hy, is the greatest field for which the curve of G 
exhibits a horizontal gradient at any point. Thus H,, is the largest field 
for which a solution is possible of the equation 


ead 
air a 
Fay reek are awe 


namo ((. 


Finally we may observe from fig. 2 that H,, is the field for which the 
curve of @ exhibits a horizontal gradient when w=0. Hence, from (1), 


Pe. fo 


It should be realized, of course, that the functions F and K in fact 
depend on temperature and on the size and shape of the particle as well as 
on w, though only w is variable in any given experiment. The relative 
values of H,, Hy, and H,, therefore depend on a variety of conditions, and 
in suitable circumstances they may coincide. No magnetic hysteresis 
would then be observed, and this is the case illustrated in fig. 4. 


and consequently 
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In order to estimate the extent of the hysteresis to be expected under any 
given conditions it is necessary to choose suitable forms for the functions 
F(w) and K(w). At present the theory of superconductivity is not 
sufficiently advanced for it to be possible to say whether: the general 
assumptions of § 4 are likely to be justified, so that there is no certainty 
that the problem can indeed be formulated in terms of these two functions. 
Still less is it possible to state the forms of F and K. However, one 
phenomenological theory has been proposed by Gorter and Casimir (1934), 
which exhibits some of the thermodynamical properties of super- 
conductors in an idealized form. In this theory, w is interpreted as the 
fraction of electrons which are in a superconducting condition, and the free 
energy per unit volume, in the absence of a field, is written in the form 


fey 
fa— yz lot2PvV(1—e)], 


where H, is the critical field at 0° K of a massive specimen, and ¢ is the 
reduced temperature, 7/7... Thus for a spherical specimen of radius a, 
F=— 8A, 2[w+ 27./(l—w)]- . . . . . (6) 
It has been found (Daunt, Miller, Pippard and Shoenberg 1948) that if the 
interpretation of w as the fraction of superconducting electrons be used 
to predict the temperature variation of the penetration depth A, by 
writing A=A)/\/w, Ay being the penetration depth at 0° x, then the 
prediction is in good accord with experimental data so far as the latter go. 
If then we assume the London penetration law the form of K(w) may be 
derived immediately. For example, for a sphere of radius a (London 1936) 


ihe 2 3 
Ka 1— sh (a4/ = oth Le 
ne ( POU coth («4/w)+ =). } (7) 


where «=a/Ao. 
Using eqns. (6) and (7) we may now give analytical form to the 


expressions (3), (4) and (5) for the characteristic field strengths, H,, 
Ay, and He rm 

fit. 2 w— 2¢[1— /(l—w )] 

H. Maxey eke Ey ie ee 

Gale: a fax 1 (3/a./w) coth (a /w)+ pu tthe; (toa ES 


)=F Max {o* 1— #/4/( (l-w 

H) oreo a4/w coth (a4/w)+a fo dosed Va) 8 3} ~v3K9) 
ebroy OAK ; 

Be) al) 0 dee eT) 


A typical example of the variation with temperature of these 
characteristic fields is shown in fig. 6, which has been computed for x 
equal to 3, i.e. a spherical particle whose radius is 3A). It will be seen 
that the hysteresis becomes less marked as the temperature is raised, 
and that for values of t greater than 0-88 the transition is reversible. The 
variation with « of the temperature, ¢’, above which no hysteresis should be . 
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observed, may be deduced explicitly from (8), since, according to § 2, it is 
only below ?¢’ that w’, the equilibrium value of w, is different from ve in a 
field H.. Hence if we put t=?’ in (8) and expand the quantity in brackets 
as a Series in w, the desired relation between ¢’ and « is obtained by equating 
to zero the coefficient of the term in w, since this ensures that the right 
hand side just achieves its maximum value when w—0. Performing 
this expansion we find that (8) takes the form 


7) — Qe? 
ee Mazes 1 (A) 
(7: a Max { am 51 Je 0 90 } , 


2:0 


a 


Magnetic field in units of H, 


Critical fields versus reduced temperature : 
(a) H, for bulk material, 
(6) Hs. for colloidal sphere of radius 3A, 
(c) H, for colloidal sphere of radius 3Ab, 
(d) Hg for colloidal sphere of radius 3A). 


Hence 
t/2 Dey? 
Sea ear ay) 
a(t?) <2] 
or 
21¢'2 
—S oe Oe e e ° e e e ll 
te 6( le) ey 


Under these conditions (H,/H,)?—=[10(1—t’?)]/«2, which: has the same form 
as (10), so that, as expected, there is no supercooling. A similar expansion 
of (9) yields the result that the superheating also just vanishes at this 


temperature. 
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Inspection of (11) shows at once that there is no size of particle, however 
small, for which some hysteresis should not be observed at a sufficiently low 
temperature. This temperature may, however, be so low as to lie outside 
the convenient experimental range. For example, the value of « for 
Shoenberg’s (1940) smallest colloid preparation was about 0-35, so that. 
from (11), we should expect t’ to be about 0-21, corresponding to a tempera- 
ture of 0-87°x. It is therefore not surprising that no hysteresis was 
observed at temperatures above 1° kK. 


§6 


The experiments of Shoenberg are the only ones so far performed on 
spherical particles, and therefore susceptible of comparison with the above 
results. It is a simple matter to modify the theory to apply to thin 
films, for which there is more information, but it is unlikely that any 
films have yet been studied which have been so uniform as to make 
detailed comparison with theory worth while, and, as has been discussed 
elsewhere (Pippard 1951), it is possible that the magnetic behaviour of 
films is explicable in terms of phase-boundary formation rather than a 
homogeneous transition of the type dealt with here. We shall therefore 
attempt a comparison of theory and experiment only for spherical particles. 

The smallest colloidal particles studied by Shoenberg had diameters less 
than A, and showed no hysteresis. This we found in the last section to be 
in accordance with the theory. It is not, however, a very convincing 
test of the theory, which can be more critically assessed on its predictions 
concerning the hysteresis of larger spheres. Unfortunately no detailed 
data are available for Shoenberg’s preparation Hg B, which showed 
greater hysteresis at lower than at higher temperatures, as the theory 
predicts, while Hg C, in which the particles were probably rather larger, 
shows the greatest hysteresis at the higher temperatures. This suggests 
that the particles are now large enough to allow w to vary from point to 
point, and their behaviour is therefore not directly comparable with our 
theory. Nevertheless a rough comparison is possible ; for the condition we 
have laid down, that w shall be constant, is an extreme condition, and the 
hysteresis observed should not exceed that predicted by theory if an 
alternative mechanism for the transition exists. 

We may estimate the mean value of « for the preparation Hg C4 from 
the fact that at the lowest temperatures the ratio y/x) of the measured 


susceptibility to that appropriate to perfect diamagnetism was 0-50. This 
means that | 


3 3 
| — — th 35 — Ure 
aoe a+ x2 0-50, 


or a=4-75, 


According to (11), then, we should expect no hysteresis for values of ¢ 
greater than 0-945, that is, at temperatures above 3-94°K. Yet experi- 
mentally it was found that at 3-883° K the hysteresis was very great, the 


ee 
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ratio H,,/H,, being nearly 3; at this temperature eqns. (9) and (10) 
predict a ratio of only 1-03. 

This striking discrepancy may reveal the inadequacy of the theoretical 
model, but another interpretation is conceivable, namely, that the 
estimate of the amount of mercury present in the preparation Hg C4 was 
seriously in error, and that in fact y/, was considerably greater than 0-50. 
There are two pieces of evidence which favour this view. First, if we 
take «—=4-65 and \),=4-3x 10-8 em, the value determined by Laurmann 
and Shoenberg (1949), the mean radius of the particles becomes 2 x 10-5cm, 
while a microscopic examination of the preparation suggested that it was 
about 10-*cm. Secondly, if the penetration depth varies with temperature 
according to the law, 

(es ele) ee ae ae) 
which describes the experimental results both for macroscopic specimens 
and for the smallest colloid preparation (see Laurmann and Shoenberg 
1949 ; Daunt, Miller, Pippard and Shoenberg 1948), the susceptibility of a 
preparation for which «4-65 should decrease by a factor 2-1 between 
0° K and 3-883° kK, while the measured decrease was only by a factor 1-24. 
This discrepancy may be accounted for if we assume that the mercury 
present in the preparation was only about 0-62 of the estimated amount, 
and that in consequence y/y, at 0° K was 0-8 rather than 0-5. The drop by 
a factor 1:24 between 0° K and 3-883°K is then consistent with (12). 
Moreover, the value of « is now raised to 15, corresponding to a mean 
radius of 6-5 x 10~-> em, which is much closer to the microscopic estimate. 
It is not immediately obvious that such a gross error in the chemical 
estimation of the mercury content of the preparation could have arisen. 
Perhaps a substantial proportion was present not in the form of metallic 
mercury but as oxide, either on the surface of the particles or as detached 
fragments, which would not become superconducting but would be 
estimated in the chemical analysis. 

If this interpretation is correct, the discrepancy between experiment and 
theory is much less, as is shown by the comparison in the table, for which 
Ay, and H,, were calculated taking «=15 and H)=410 gauss. 


IMA H sh ta § sh jek sc Jel sc 

cale. obs. cale. obs. 
1-935 690 410 Wal 210 
3°883 75 7(5) 31 28 


The observed values of H,, and H,, have been taken as those field strengths 
at which about one half of the material had made the transition. ‘The good 
agreement between theory and experiment at the higher temperature is 
probably largely fortuitous, in view of the doubtful interpretation of the 
experimental results. At the lower temperature agreement is not nearly so 
good, but at least the disagreement is in the right direction, the observed 
hysteresis being less than the predicted. It is quite possible that the 
hypothesis made in § 2, that w is constant over the particle, begins, as the 
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particle size is increased, to break down at the lower temperatures first, 
since the observed temperature variation of the interphase surface energy 
indicates that the range of order is smaller at lower temperatures. Thus 
one might expect the comparison of theory and experiment to give results 
similar to those in the table. 

It would be unwise, however, on the basis of so little experimental 
material to claim more than qualitative validity for the theory. A much 
better check could be obtained from experiments with particles for which « 
was about 3, so that hysteresis should, according to (11), be observed at the 
lower but not at the higher temperatures, as illustrated in fig. 6. Such 
experiments, if they gave results in conflict with the theory, might serve 
a useful purpose in suggesting a possible form for the dependence of the 
free energy on the degree of order, w, if indeed the thermodynamics of a 
superconductor may be formulated in so simple a way. 
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ABSTRACT 


The angular distribution of high energy neutrons produced by bombard- 
ing a beryllium target with the 170 Mey proton beam of the Harwell 
cyclotron has been measured. The angular spread for neutrons with 
energies greater than 50 Mev was found to be 67° at half intensity and in 
this spread the higher energy neutrons were shown to be confined to the 
smaller angles. The total cross section for the production of neutrons 
with energies greater than 50 Mev was found to be 99-+40 x 10-2? cm?2, 
while the differential cross section at 4° was 69 +30 10-27 em?/steradian. 


§1. LyrropucTION 


THE angular distribution of high energy neutrons produced by 340 Mev 
protons has been measured for angles up to 27° to the forward direction 
and for various target materials by Miller, Sewell and Wright (1951). 
In this work the relative numbers of neutrons were indicated by the 
UC positron activity induced in polythene discs exposed to the neutron 
beam. De Juren (1951) made similar measurements by recording fissions 
in bismuth. The excitation functions for these reactions are not exactly 
known but they certainly differ in shape; moreover, the threshold 
energy for the carbon reaction is about 20 Mev whereas the fission of 
bismuth is not detectable below 50 Mev. Consequently the very similar 
results obtained in the two experiments indicate that at 330 Mev’ the 
angular distribution is not strongly dependent upon neutron energy. 

It was thought worthwhile to extend this work using the 170 Mev 
protons accelerated by the Harwell cyclotron. A beryllium target was 
used and the neutrons produced were detected with a recoil proton 
telescope (Hadley et al. 1949). An absorber was used to vary the threshold 
of detection and the response as a function of neutron energy calculated 
with the help of the Berkeley experimental data on n-p scattering 
(Christian and Hart 1950). Owing to the complication of this method 
of measurement it was only convenient to make measurements at a few 
different angles. 

During this work the opportunity was taken to redetermine the cross 
section for neutron production in beryllium and the result is included. 


* Communicated by the Author. 
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§2. APPARATUS 


The experimental arrangement of fig. 1 was used. There were six 
collimating holes in the main 6 ft. thick concrete shielding wall W. The 
axis of hole No. 6 was tangential to the orbit of radius 48-5 in., but the | 
axes of holes 1-5 intersected such tangents at angles of 4°, 26°, 40°, 51° 
and 60° respectively. A monitor target D was aligned for hole No. 6 
and the neutron intensity in the forward direction from this target was 
measured in the double coincidence telescope M. The beryllium target 
B, 0-352 in. thick, was aligned with one of the holes 1-5 while the other four 
were blocked with brass absorbers. The neutrons which emerged at 
this angle to the forward direction passed through a thin aluminium window 


Fig. 1 


CYCLOTRON 


V in the cyclotron tank and: were then recorded in the triple coincidence 
telescope 7’. The front collimating block F' was also aligned with the 
particular collimating hole in use. 

In addition a measure of the integrated neutron output from the 
cyclotron during any particular run was obtained from the counts 
recorded by a slow neutron counter placed in a recess in the shielding 
wall at S. 

The neutrons emerging from the 5-3cm diameter collimator hole were 
allowed to strike a polythene disc R, 9mm thick and 5 ecm diameter, 
which had been set up on the axis of the collimator. Three scintillation 
counters were used to detect the protons ejected from the polythene disc 
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at an angle of 12° to the forward direction. Counters 1 and 2 were 
identical and used clear naphthalene crystals 3 cm square and 3-2 mm 
thick, whereas counter 3 used a larger crystal 4-8 em diameter and 2-5 em 
thick. A carbon absorber A could be placed between counters 2 and 3 
to set the minimum energy of the protons detected. Taylor’s (1950) 
results for the range energy relation in carbon were used. The solid 
angle subtended at the radiator by the telescope was defined by the 
2-45 cm hole in a 4 cm long copper block C placed immediately after 
counter 2. Only protons traversing this hole reached the absorber A 
and counter 3. The monitor channel was similar but used only two 
scintillation counters and was set for a threshold proton energy of 90 Mey. 


§ 3. EXPERIMENTAL MEeTrHop 
Tests were first performed to determine that every proton traversing 
the three counters of the measuring telescope was registered as a triple 
coincidence. Similar tests were made in the monitor channel. The 
measuring telescope was then set to a particular threshold energy and the 
counting rates obtained at each hole collated with the help of the monitor 
channel. This experiment was repeated for several threshold energies. 
In practice some of the recoil protons failed to reach the third counter 
due either to nuclear absorption or scattering in the absorber. The 
correction for this effect was determined by an experiment in the external 
proton beam of the cyclotron. It was found that the correction for a 
particular absorber depended on the incident proton energy. To apply 
the correction, therefore, a rough knowledge of the incident proton 
spectrum was required and this was obtained from the data already 
recorded for each neutron angle. The corrections varied from 10%, for 
the lowest threshold energy to 40% for the highest threshold energy. 
The Berkeley results for n—p scattering were then applied to the corrected 

proton data to give the incident neutron intensities. 


§ 4. RESULTS 

In fig. 2 the angular distribution of neutrons is drawn for several 
different values of the threshold neutron energy. The curves are 
normalized to unity for the forward direction and the angles in the laboratory 
system are plotted as abscissae. For neutrons of energy greater than 
50 mev the full angular spread at half intensity is 67°. The higher 
energy neutrons are confined to the smaller angles, which is the result to 
be expected if they are produced by interaction of the incident proton 
with a single nucleon within the nucleus. 

The curve for the lowest threshold neutron energy was integrated and 
the effective proton current estimated from the results of Cassels, Dickson 
and Howlett (1951). Hence the total cross section for the production in 
beryllium of neutrons with energies greater than 50 Mev was found to be 
99-40 10-27 cem?. ‘The differential cross section at 4° was found to be 
69-430 x 10-2? em2/steradian in reasonable agreement with the value of 
99 +19 x 10-2? cm2/steradian obtained earlier by Cassels, Randle, Picka- 
vance, and Taylor (1951). 
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ABSTRACT 


A preliminary survey is first presented of the factors which have been 
suggested to account for the generally observed decrease in sensitivity of 
a photographic emulsion on lowering. its temperature to that of liquid 
nitrogen. It is concluded that the main problem to be resolved is whether, 
after due allowance for the decreased absorption of the halide, the observed 
speed loss can be completely assigned to the effective high intensity 
reciprocity failure of the low temperature exposure, or whether some 
additional mechanism of exposure wastage comes into play. Experiments 
are then described, the results of which show that for a range of emulsions 
the observed speed loss on passing from room to liquid nitrogen temperature 
can be quantitatively accounted for in terms of decreased absorption of 
the halide and the fact that low temperature exposures are essentially of 
very high intensity. The fluorescence of bromo-iodide emulsions whilst 
undergoing a low temperature latent-image exposure is apparently not 
connected with the speed loss, but is merely a visual manifestation of 
the recombination which is occurring in an unobservable manner at room 
temperature. Finally the results of other investigators are discussed in 
the light of the new findings; several pieces of evidence emerge which 
confirm their generality. ' 


§1. INTRODUCTION 


THE results of experiments on the influence of temperature on photographic 
phenomena have provided confirmation of the validity of the fundamental 
hypotheses of latent-image theory now generally accepted. These were 
first clearly enunciated by Gurney and Mott, and state that latent-image 
formation takes place in two steps, namely, the release and subsequent 
trapping of an electron from a halide ion of the photographic grain, 
followed by neutralization by a silver ion, the process being repeated 
several times in order to form a silver speck capable of promoting 
development (Gurney and Mott 1938). 

It is not difficult to explain how a latent image can be formed as the 
result of an exposure at low temperatures (e.g. liquid nitrogen, 77° K), 
but considerable uncertainty has persisted as to the factors which come 
into play in order that the sensitivity of a non-dye sensitized emulsion 
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be reduced to about a tenth of its room temperature sensitivity (Sheppard, 
Wightman and Quirk 1934, Webb and Evans 1938, Berg 1939). An 
investigation of these factors was the purpose of the work now being 
reported. 

Anticipating the conclusions to be drawn from the experimental 
results it has been found that for all emulsions examined the fall in 
sensitivity on proceeding from room to liquid nitrogen temperature can 
be accounted for in terms of the decrease in absorption coefficient of the 
silver halide, and the fact that, from the point of view of the behaviour 
of the emulsion, low temperature exposures are effectively of very high 
intensity. In order to bring out the precise aim of the experiments 
performed and the meaning of the results, it is first necessary to review 
briefly previous suggestions as to why the sensitivity of a photographic 
emulsion is in general observed to decrease on temperature reduction. 


§2. PRELIMINARY SURVEY 


The most obvious effect of temperature variation is an alteration in 
the absorption coefficient of the silver halide of an emulsion. Toy 
and Harrison have shown that the absorption coefficient of pure 
silver bromide for the wavelengths of 43504, 40474 and 36504 falls 
by a factor of three on reduction of temperature from 20° c to —180°¢ 
(Toy and Harrison 1930). From the results of direct absorption 
measurements on photographic emulsions by Biltz it can be deduced that 
the speed of a pure bromide emulsion to tungsten radiation should be 
reduced by about 0-5 log H units on lowering the temperature from room 
to liquid nitrogen (Biltz 1949). For a bromo-iodide emulsion containing 
about 5 mol. per cent of silver iodide the corresponding speed loss should 
be about 0-4 log # units (Biltz 1951). 

A second suggestion, stimulated by the nature of experimental results 
pertaining to the phenomenon of high intensity reciprocity failure, was 
that low temperature exposures were essentially of very high intensity. 
Berg has shown that at room temperature the reciprocity failure curve 
of an emulsion becomes flat below a time of 4 10~-° seconds, i.e. below 
this time the emulsion demonstrates a constant high intensity failure 
speed loss (Berg 1940). On lowering the temperature of the emulsion 
this flat portion of the reciprocity failure curve commences at longer 
exposure time. At liquid nitrogen temperature emulsions do not show 
any reciprocity failure (Berg 1939), and it is reasonable to suppose that 
the observed flat reciprocity failure ‘curve’ is merely the flat high 
intensity portion observable at room temperatures shifted to cover the 
whole range of normal exposure times. If this view is correct the speed 
loss observed on passing from room to liquid nitrogen temperature, if 
measured from the optimum of the room temperature reciprocity failure 
curve, will contain a component representing the speed loss which ensues 
on passing from the optimum to the exposure level at which the reciprocity 
failure curve flattens out, and which should in theory be capable of 
measurement at room temperature. 
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The general opinion to be derived from the literature on the subject 
of low temperature sensitivity is that decrease of absorption coefficient 
and effective high intensity reciprocity failure are insufficient to account 
for the observed loss of speed on temperature reduction (Berg 1948). 
It has always. been assumed that some other mechanism of exposure 
wastage comes into play. The observation by Meidinger that bromo- 
iodide emulsions fluoresced strongly at liquid nitrogen temperature 
appeared to provide an explanation of the extra speed loss (Meidinger 
1939). It was supposed that the number of electrons lost by recombination 
at low temperatures was greater than the number lost at room temperature, 
this loss being manifested by the fluorescence of the halide. Alternatively 
it could be said that no electrons were lost in a grain at room temperature, 
while some were lost at low temperatures. It has been observed by 
Farnell and Burten (1951) that a bromo-iodide emulsion fluoresces with 


Fig. 1 
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Illustration of the interrelation between room- and low-temperature sensitivity 

; in terms of a reciprocity-failure diagram. BK is the room temperature 
reciprocity-failure curve. The total speed loss, S, is equal to the total 
high-intensity failure, H, plus the quantity, A, which may represent 
reduced absorption alone, or may include a factor covering increased 
recombination. 


the characteristic green colour of silver bromo-iodide (Meidinger 1939, 
Farnell, Burton and Hallama 1950) during a latent-image exposure to 
‘ultra-violet radiation. On the other hand it was not proved that this 
observed recombination was either merely a visual manifestation of that 
occurring at ordinary temperatures, or that it represented additional 
recombination over that occurring at room temperature. 

The general situation regarding the loss of speed observed on 
temperature reduction as outlined in this survey may be conveniently 
summarized in diagrammatic form as illustrated in fig. 1. ABR is the 
reciprocity failure curve of an emulsion at ordinary temperatures, the 
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point B representing the ‘bend over’ point first observed by Berg. 
LL is the ‘ reciprocity failure curve ’ of the emulsion at low temperatures, 
which may be regarded as a portion of RB shifted downwards by decrease 
in absorption, and possibly by increased recombination, and to longer 
times. The total log speed loss (S) from the room temperature optimum, 
is thus composed of the high intensity failure portion (), plus a portion A 
which may be due to a fall in absorption alone, or may in addition contain 
a factor representing increased recombination at the low temperature. 


§3. PRINCIPLE OF THE EXPERIMENTS 


It will be appreciated from the previous discussion that the problem to 
be resolved in the field of low temperature sensitivity is whether, after 
allowance has been made for decreased absorption, the effective high 
intensity reciprocity failure of a low temperature.exposure is sufficient 
to account for the observed speed loss, or whether other forms of exposure 
wastage come into play, for instance, increased recombination. The 
obvious experimental approach is to determine the total speed loss 
from the room-temperature optimum to a low-temperature exposure, and 
compare it with the sum of the total speed loss due to high-intensity 
failure (from the optimum to the bend-over point) and the speed loss due 
to fall in absorption. If the total loss exceeds the sum of the latter 
effects then it is reasonable to suppose that some additional mechanism 
of exposure wastage is coming into play. It is to be noted that this 
approach is determined by experimental facts, and the results will be 
valid whatever the explanation of high-intensity reciprocity failure and 
the bend-over point. 

In practice several deviations from the above scheme were made in 
order to overcome experimental difficulties. For instance, it is impossible 
in the present state of our knowledge to calculate precisely the speed 
variation due to change in absorption of an emulsion for a band of 
wavelengths in the visible region, as is employed in exposures to tungsten 
radiation. This difficulty was largely resolved by making exposures to 
mercury 3650 A radiation, which was absorbed to the extent of about 99% 
at ordinary temperature and 95% at low temperatures, with the result 
that observed speed losses on temperature reduction were essentially 
due to causes other than change of absorption. Some quantitative data 
on this topic will be presented later. The smaller dependence of sensitivity 
on temperature at short wavelengths has been demonstrated by Webb 
(Webb 1935). 

A second difficulty was that the shortest exposure time possible with 
the available reciprocity sensitometer was only 1-58 10-4 sec, which 
is more than ten times longer than the bend-over time at ordinary 
temperatures. But since the high-intensity portion of the reciprocity 
failure curve is shifted approximately 0-6 logt units towards longer 
times for a decrease in temperature of 10°, it was a simple matter to 
record exposures at and beyond the bend-over point when the emulsion 
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was at the temperature of solid carbon dioxide (—78-5° c). Details of 
the experimental technique and emulsions employed will now be 
presented. 


§4. EXPERIMENTAL DETaAILs AND TECHNIQUES 
4.1. Apparatus and Exposure Technique 


A reciprocity sensitometer, with the normally fitted plateholder removed, 
was employed to give step wedge exposures to 6 in. <1 in. strips of film 
clipped on the inside wall of a Pyrex Dewar flask. Comparative exposures 
of a particular emulsion were made with the flask empty and full of 
liquid nitrogen, the exposure at room temperature being made at the 
optimum of the previously determined reciprocity failure curve. It was 
sometimes impossible to do this with the fairly low intensity of ultra- 
violet radiation available, and so the results were adjusted to the optimum 
from the reciprocity-failure curve. The strips of film were backed by a 
stout strip of plastic to prevent them buckling at the low temperature. 

Extraneous radiation was prevented from reaching a strip undergoing 
exposure by sheathing the Dewar flask in black paper, except for an 
exposure slot of dimensions slightly larger than the step wedge. 

Reciprocity failure curves at the temperature of solid carbon dioxide 
(—78-5° Cc) were obtained by exposing strips in a similar manner to that 
used at liquid nitrogen temperature. Cooling was achieved by immersing 
strips in the Dewar flask full of ethyl or methyl alcohol with the lower 
half occupied by lumps of solid carbon dioxide. It was proved by several 
control experiments that immersion in alcohol had no influence whatsoever 
on the sensitivity of the emulsions employed. 

At least two strips, and quite often four, were exposed in order to 
obtain any particular characteristic curve ; in all cases good reproducibility 
was obtained. The ordinary lamps of the sensitometer were employed for 
tungsten radiation exposures. A mercury 125 w high-pressure are with 
a Wratten 18A filter was employed to give ultra-violet exposures 
predominantly of wavelength 3650 a. 


4.2. Emulsions 

In order that emulsions could be examined at different stages of 
sensitivity, as determined by their preparation, and free from impurities 
normally present in commercial emulsions, work was confined. to simple 
experimental emulsions. The basic halide compositions were pure silver 
bromide and silver bromide-+_5 mol. per cent silver iodide, ripening being 
carried out in excess halide solution. Details of emulsion preparation, 
together with the reference letters employed in this paper, are given in 
table 1. None of the emulsions were colour sensitized. 


4.3. Development 

Since the aim of this work was to achieve some form of general result 
no particular attention was paid to the type of development employed. 
Bromo-iodide emulsion strips were developed. for 16 min in Kodak D19b 
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developer (and occasionally for shorter times), and pure bromide strips 
for 5 min or 6 min. All the emulsions, with the exception of I (sulphur 
sensitized) formed predominantly internal image. Because of this fact 
no complications due to solarization appear to have been encountered. 


§5. Resuxtts, DiscussIon AND CONCLUSIONS 


' 


5.1. Absorption Change on Temperature Reduction 


The aim in this section is to show that from the recordings of room- 
and low-temperature characteristic curves to tungsten and ultra-violet 
radiation a reasonable value for the speed loss to tungsten radiation 


Table 1. Treatment of Experimental Emulsions 


AgBr AgBr-+5 mol. per cent AgI 


A. Ripened 15 min at 40° c ; . Ripened 5 min at 40° c; 
undigested. undigested. 


B. Ripened 45 min at 40° c ; . Ripened 15 min at 40° c ; 
undigested. undigested. 


C. Ripened 45 min at 40° c ; . Ripened 45 min at 40° c ; 
digested 55 min at 55° c. undigested. 


G. Ripened 15 min at 40° c ; 
digested 35 min at 55° c. 


H. Ripened 45 min at 40° c ; 
digested 35 min at 55° c. 


. Ripened 15 min at 40° ¢ ; 
digested 35 min at 55° c, 
in presence of a sulphur 
sensitizer. 


due to absorption change may be detected. Some typical results for a 
bromide and a bromo-iodide emulsion are illustrated in figs. 2 and 3 
respectively. Section (a) of each figure compares the room- and low- 
temperature characteristic curves to tungsten radiation, while section (b) 
effects the same comparison for ultra-violet radiation. It will be noted 
immediately that the speed losses for ultra-violet radiation are considerably 
less than those for tungsten radiation, this being attributed to the fact 
that the fraction of incident tungsten radiation absorbed falls considerably 
on temperature reduction, whereas there is little change in the fraction 
of ultra-violet radiation absorbed. In fig. 4 a further comparison for 
the bromo-iodide emulsion is made; in this case a short development 


———— 
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Fig. 2 


DENSITY 


DENSITY 


05 


20 


() 

Room- and low-temperature characteristic curves (marked F# and L respectively) 
for a pure bromide emulsion (C). (a) to tungsten radiation, (6) to ultra- 
violet radiation. The pure bromide emulsion’was developed for 6 min in 
D19b. The densities have been corrected for fog. 
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Fig. 3 
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1-0 DENSITY 


(0) 
Room- and low-temperature characteristic curves (marked R and L respectively) 
for a bromo-iodide emulsion (I). (a) to tungsten radiation, (b) to ultra- 


violet radiation. The*bromo-iodide emulsion was developed for 16 min in 
D19b. The densities have been corrected for fog. 
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Fig. 4 
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(0) 
Room- and low-temperature characteristic curves for the same bromo-iodide 
emulsion (I) on short development (4 min D19b), (a) to tungsten radiation, 
(6) to ultra-violet radiation. 
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was employed. The actual speed losses are much larger than for the 
long development, but the tungsten speed loss is still greater than the 
ultra-violet speed loss. This result demonstrates that speed loss on 
temperature reduction involves a process of the high-intensity failure 
type, since prolongation of development is a recognized means of reducing 
high-intensity failure. 

If it is assumed that absorption changes have little influence on the 
speed of an emulsion to ultra-violet radiation, then it will be clear that 
the difference of the speed losses to tungsten and ultra-violet radiation 
respectively should be equal to the speed loss due to reduced absorption 
of tungsten radiation at the low temperature. Since the speed losses 
are measured from the optimum of the room-temperature reciprocity 
failure curve, high-intensity failure or other forms of speed loss will be 
the same for both tungsten and ultra-violet exposures. A constant 
difference between tungsten and ultra-violet speed losses is to be expected 
if the emulsion coating weights are constant. This latter point was 
checked by silver determinations and by measuring the optical densities 
of unexposed strips; the maximum variations in coating weight were 
found to be about 20%. In table 2 the speed losses to tungsten and 
‘ultra-violet radiation are quoted for densities of 0-2 and 0-7 above fog. 


Table 2. Speed Losses on Temperature Reduction from 20° c to —196° 6, 
: for Reference Densities of 0-2 and 0-7 above Fog 


Values quoted in log # units 


D=0:2 above fog D=0-7 above fog 
Emulsion and time of ©|_————@—@————_—_____+—_|_ 
development in DI9b |Tungsten| U.V. Tungsten} U.V. 
speed |speed| Difference] speed |speed| Difference 
loss | loss loss loss 
Bromide A. 5 min 0:78 | 0:20 0-58 0-85 | 0-14 0-71 
Bromide B. 6 min 0:89 |0:12| . 0-77 0-91 0:07 0-84 
Bromide C. 6 min 1:07 0-34 0-73 0-75 0-18 0-57 
Bromo-iodide D. 5 min 1:95 | 1-42 0°53 — — = 
Bromo-iodide D.16 min] 1-23 | 0-90 0:37 1:33 | 0:83 0-50 
Bromo-iodide E. 16 min 1-10 | 0-68 0-42 0-94 | 0-51 0-43 
Bromo-iodide F. 16 min 1:23 | 0:79 0-44. 0:99 0-64 0-35 
Bromo-iodide G. 16 min 1-34. 0-86 0:48 1-41 0-74 0-67 
Bromo-iodide H. 16 min 1-07 0-76 0-31 0:89 | 0-53 0-36 
Bromo-iodide I. 4 min 1-53 123 0-30 1-78 1:34 0-44 
Bromo-iodide I. 16 min 0-71 0-40 0-31 0-79 0-34 0:45 


It will be noted immediately that the difference factors for pure bromide 
emulsions are larger than those for bromo-iodide emulsions. This is 
consistent with the results of Biltz mentioned earlier, namely, that the 
decrease in absorption of a bromo-iodide emulsion (~5 mol. per cent AgT) 
on temperature reduction is somewhat less than that of pure bromide. 
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The difference values for the two different types of emulsion are as 
constant as can reasonably be expected in view of the variations in 
coating weight, and the experimental errors which accrue in the 
derivation of the speed loss difference. . 

In the case of the bromo-iodide emulsions, for which more extensive 
results are available, the speed loss difference appears less, on the average, 
at the lower reference density. 

This variation with reference density can be understood by referring 
to an expression for the relative log speed of an emulsion derived by 
Burton (unpublished work), and confirmed by checking against experi- 
mental data (unpublished work by Farnell). The expression is : 


, log;9 S=log,, d—xA-+const., . . . 2... (1) 
where S=speed in reciprocal quantic units, 4—density of unexposed 
emulsion to the wavelength of the exposing radiation, and x=ratio of 
reference density to maximum density. From expression (1) a relation 


for the change in log speed on change of emulsion density can be derived, 
and is quoted below, . 


10g sp S100 15.(43) 45) e(A vd.) te oe (2) 


Ifthe change in emulsion density which occurs on lowering the temperature 
is such that at a low reference density the first term on the right hand 
side exceeds the second, raising the reference density (a) will diminish 
the speed change which occurs on lowering the temperature. Thus 
it appears that the smaller speed loss difference (tungsten—ultra-violet) 
at the lower reference density can be assigned to a small absorption 
speed loss with the ultra-violet exposures, which effectively reduces the | 
difference between the tungsten and ultra-violet speed losses. 

This deduction was checked by measuring the density of unexposed 
emulsions to 3650 A radiation at room- and low-temperatures and inserting 
the values in expression (2). The densities were determined by exposing 
two strips one behind the other, developing them together, and measuring 
the separation of the characteristic curves along the log exposure axis. 
In this way it was found that the densities of the bromo-iodide emulsions 
at room temperature were on the average 2:3 and fell to 1-5 at liquid 
nitrogen temperature. The average maximum density was about 3-0. 
Insertion of these values in expression (2) shows that for 3650 a radiation, 
an absorption speed loss of 0-16 log # units occurs at a reference density 
of 0-2 above fog, but. that no absorption speed loss occurs at a reference 
density of 0-7 above fog. The quantitative agreement of this result 
with the discrepancy pointed out above was taken as confirmation of 
the explanation, and also as a confirmation of the original assumption 
that speed losses due to change in absorption of ultra-violet radiation 
are small. 

To sum up this section, it may be said that comparisons of the speed 
loss to tungsten and ultra-violet radiation under equivalent conditions 
appear to provide a satisfactory method of determining the absorption 
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speed loss with the former type of radiation. For an average bromo-iodide 
emulsion this quantity is about 0-4 log Z units, and about 0-6 log # units 
for a pure bromide emulsion. 


5.2. High Intensity Reciprocity Failure 


It has now been shown that comparisons of room- and low-temperature 
sensitivity to ultra-violet radiation are sensibly free of any complication 


due to absorption change at a density of 0-7 above fog; a reasonable — 


estimate can be given of the absorption speed loss at 0-2 above fog. In 
order to complete the correlation between room- and low-temperature 
sensitivity it remains to be determined whether the observed ultra-violet 
speed losses can be completely assigned to high intensity reciprocity 


failure of the low-temperature exposure, or whether there is some_ 


additional mechanism of exposure wastage coming into play. 

Typical reciprocity failure curves for a bromide and bromo-iodide 
emulsion at room temperature (20°c) and ¢@solid carbon dioxide 
temperature (—78-5° c) are shown in figs. 5 and 6 respectively. The 
emulsions are actually those for which the characteristic curves have 
been presented (silver bromide and silver bromo-iodide). As described 
earlier the —78-5° © reciprocity failure curve is the very high intensity 
portion of the room temperature curve shifted to longer times; the 
‘ bend-over ’ points are well defined. Owing to the rather low sensitivity 
of the experimental emulsions employed in this work, and the large drop 
in sensitivity which occurred on reduction of the emulsion temperature 
to —78-5° c, it was only possible to obtain satisfactory reciprocity failure 
curves at this temperature for densities of 0-2 above fog and with long 
developments. An exception, however, was the more sensitive sulphur 
sensitized emulsion for which satisfactory reciprocity failure curves were 
obtained at 0-7 above fog. 

In table 3 comparisons are made of the ultra-violet speed losses on 
passing from 20°c to —196° c, with the high-intensity failure speed 
loss reckoned from the room temperature optimum to the horizontal 
portion of the —78-5°c curve. The latter curves were shifted up 
0-12 log # units to counteract the small reflection loss which occurred 
at the exposure slot in making the exposures; the room temperature 
curves were obtained with the step wedge in direct contact with the 
sheet of film undergoing exposure. 

A good correlation between the values in the second and third columns 
can be observed. The ultra-violet speed losses are invariably slightly 
less than ‘the high-intensity failure speed losses. The reason for this 
becomes apparent when it is realized that the former values are over- 
estimated, for the purposes of this comparison, by an amount corres- 
ponding to the speed loss due to the fall in absorption of ultra-violet 
radiation on temperature reduction from 20° ¢ to —196°c. The latter 
values are over-estimated by a larger amount corresponding to the speed 
loss due to the fall in absorption of tungsten radiation on proceeding 
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Pe a anes The bracketed values in the second column are 
i speed. losses at 0-2 above fog, corrected for absorption 
ic kale eee emulsion densities to 36504 radiation at room- 
Fe. peratures and applying expression (2). From a consideration 

e results detailed in table 2 and of the temperature ranges involved 
the high-intensity failure speed losses noted in the third column of 


Figs. 5 and 6 


2-6 


2:8 


3-4 


3-6 


4 @ > ! @) | 2 


Reciprocity-failure curves of the bromide emulsion (C) and the bromo-iodide 
emulsion (I) at 20° c and —78-5° c (5 and 6 respectively). The reference 
density of the bromide emulsion curves is 0-2 above fog ; for the bromo- 
iodide emulsion, curves are given for both 0-2 and 0-7 above fog. 
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table 3 should exceed the corrected ultra-violet speed losses by about 
0-2 log B units for bromo-iodide and 0-3 log H units for pure bromide 
emulsions. These values are the estimated absorption corrections for 
tungsten radiation in passing from 20°c to —78-5° co, the reasonable 
assumption being made that the speed loss is approximately linearly 
related to the absorption change. The above condition will be seen to 
have been fulfilled within the limits of experimental error. 


Table 3. Comparisons of Ultra-violet and High Intensity Reciprocity 
Failure Speed Losses 
Reference density 0-2 above fog. All values in log # units 


Emulsion and development Ultra-violet High intensity 
speed loss failure speed loss ° 
Bromide B. 6 min D19b 0-12 (—0-09) 0-21 
Bromide C. 6 min D19b 0-34 (+0-17) 0-42 
Bromo-iodide E.16 min D19b | 0-68 (0-54) 0-76 
Bromo-iodide F.16 min D19b | 0-79 (0-67) 0-87 
Bromo-iodide G.16 min D19b | 0-86 (0-74) 0-90 
Bromo-iodide H. 16 min D19b | 0-76 (0:58) 0-79 
Bromo-iodide I.16 min D19b | 0-40 (0-28 0-41 : 
; 0-34 (0-7 above fog) | 0-50 (0-7 above fog) 


It appears, therefore, that the observed ultra-violet speed losses can 


be completely accounted for in terms of the high-intensity reciprocity. 


failure of the low-temperature exposures. The generality of the result is 
confirmed by the fact that the pure bromide emulsions suffered little 
speed loss on temperature reduction, and show little high-intensity 
failure, while the sensitivity of the bromo-iodide emulsions fell con- 
siderably and was paralleled by strong high-intensity failure. Some 
check experiments confirmed the fact that emulsions D and I showed 
very strong high-intensity failure on short development. Although the 
bend-over point could not be reached at —78-5° c, the total high-intensity 
failure visible was almost sufficient to account for the much larger speed 
losses to ultra-violet radiation on short development. 
5.3. The Fluorescence of Photographic Emulsions 

The method of exposing strips at low temperatures enabled them to 
be observed during the actual latent-image exposure. When exposed to 
ultra-violet radiation pure bromide emulsions fluoresced a faint blue colour 
due to their gelatin content, but bromo-iodide emulsions invariably 
fluoresced with the characteristic green colour of bromo-iodide mixed 
crystals. ¢ 

The interesting point about this fluorescence during exposure is its 
strong indication that the quantum efficiency of formation of latent-image 
silver is less than unity. The argument is as follows. It has been shown 
that the speed loss of a photographic emulsion on temperature reduction 
is due to reduction in absorption and high-intensity reciprocity failure ; 
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no additional mechanism appears to come into play. Now it is generally 
considered that inefficiency at high intensities is merely due to a greater 
dispersion of the latent-image silver, and does not involve increased 
recombination of electrons and positive holes. Confirmation of this 
point is provided by the fact that reciprocity failure disappears at very 
high intensities. Thus the recombination which is observed to take place 
at low temperatures must also occur at room temperature, and, as such, 
represents a quantum. efficiency of silver formation of less than unity. 
Actual measurements of quantum efficiency which involve the production 
of much larger quantities of silver, for instance the recent experiments 
of Meidinger, have provided results which lead to a similar conclusion 
(Meidinger 1949). 

5.4. Results of Other Investigators and Conclusions 

Since the results described appear to provide a general explanation of 
the normally observed loss of sensitivity on temperature reduction, it is 
of interest and importance to determine whether the findings of other 
workers are consistent with the present viewpoint. First it will be observed 
that the normally reported speed loss of ~1-0 log H units between the 
room-temperature optimum and the low-temperature exposure can be 
easily provided by the sum of the absorption and_ high-intensity 
reciprocity failure effects. Meidinger has shown that no obvious correlation 
exists between loss in sensitivity on temperature reduction and the 
fluorescence of the halide at the low temperature (Meidinger 1940). The 
results described in this paper confirm and extend his observations by 
showing that though low-temperature exposures are partially wasted by 
recombination, there is no necessity to invoke this effect in order to 
account for the observed loss’ in sensitivity. In fact it appears that 
recombination ‘takes place to equal extents at both room- and low- 
temperatures. : 

The loss in sensitivity observed on proceeding from liquid-air to liquid- 
hydrogen temperature (90° K to 20°K) is very much smaller than the 
corresponding quantity observed on passing from room- to liquid-air 
_ temperature (Berg and Mendelssohn 1938). Since reciprocity failure is 
absent at the two low temperatures the loss in sensitivity can only be 
assigned to reduced absorption. From the data already presented it 
can be deduced that for an emulsion of normal thickness the speed loss 
on reducing the temperature from that of liquid air to that of liquid 
hydrogen should be approximately 0-2 log # units. Characteristic curves 
of Process film recorded by Berg and Mendelssohn at the two above- 
mentioned temperatures show that there is indeed a loss in sensitivity 
of about 0-2 log # units, at least for low and medium densities. Thus 
it appears that even at liquid-hydrogen temperature no other factors 
leading to lower efficiency are encountered other than reduced absorption 
and high-intensity reciprocity failure. . 

The large drop in sensitivity of the surface image of chemically 
sensitized emulsions on reduction of temperature, as reported by Biltz 
(Biltz 1949) is readily explained as being due to the well-known strong 
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high-intensity failure of the surface image in such emulsions. It has 
already been described in § 5.2, how large speed losses observed on short, 
predominantly surface, development are very much . reduced when 
development is prolonged. Biltz explained his results by saying that 
chemical sensitization becomes ineffective at low temperatures. This 
form of words may be accepted if it is added that the same breakdown 
occurs at room temperature and short exposure times, but the phrase 
is unfortunate since the speed of the chemically sensitized emulsion at 
low temperature did not fall simply to the speed level of the corresponding 
unsensitized emulsion, but considerably below it. 

A similar viewpoint is applicable to the experimental results of Webb 
and Evans that the gamma of a low-temperature characteristic curve 
approaches that of an ordinary temperature curve as the number of 
low-temperature exposure intervals followed by warming up in the dark 
was increased, the total exposure in each case being kept constant (Webb 
and Evans 1938). A residual loss of speed was observed, however, which 
is to be assigned to reduced absorption. The practice of warming up 
after a section of the low-temperature exposure, and thus allowing some 
silver to be formed on which electrons can be trapped in the next section 
of the exposure, is just that required to’‘eliminate high-intensity reciprocity 
failure ; the latent-image silver is not deposited in a disperse manner as 
after a single low-temperature exposure, but is allowed to build up as a 
few large specks. Arguing on the basis of the results described in this 
paper it seems certain that a result identical with that of Webb and 
Evans could have been obtained by giving strips single low-temperature 
exposures, and processing them for increasing lengths of time in a 
commercial developer. The mechanism coming into play in this case 
would be different to that involved in the experiment of Webb and 
Evans, but the final result, viz., the reduction of high-intensity reciprocity 
failure, would be the same. In order to avoid possible confusion it should 
be pointed out here that high-intensity failure invariably contains a 
component representing loss of gamma, and in some cases represents 
this effect alone. 

A further observation now readily explained is that the sensitivity of 
an emulsion to X-rays is little temperature-dependent (Eggert and 
Kleinschrod 1941). The transit time of the photoelectron through a 
grain is so short that its influence corresponds to that of an extremely 


brief light exposure. Thus in so far as an emulsion exposed in the above - 


manner is demonstrating maximum high-intensity failure at room 
temperatures, reduction of temperature can have only small effects 
related to the slightly increased energy required to free bound electrons 
in the silver halide. . 

Finally it is worth describing some experimental results of Berg, 
Marriage and Stevens which provide excellent confirmation of the 
deductions made from results described in this paper (Berg, Marriage 
and Stevens 1941). They gave spark exposures to strips of a medium- 
speed blue-sensitive emulsion at both room and low temperatures and 
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processed strips in surface and internal developers (the latter after 
bleaching of the surface image). The surface image was considerably 
less sensitive than the internal but the speed loss on temperature reduction 
was approximately the same for the two types of image. This is significant 
because the spark exposure is of such short duration that it corresponds 
to an exposure beyond the bend-over point, and thus a decrease of 
temperature should only decrease the sensitivity by an amount corres- 
ponding to the decreased absorption. In turn the latter quantity should 
be independent of the method of development employed, which is indeed 
experimentally observed. 

Thus in conclusion it can be said that as regards the straightforward 
question of the interrelation of room- and low-temperature sensitivity of 
non-colour-sensitized emulsions, all the available evidence is consistent 
with the view that the normally observed fall in sensitivity is due to the 
high-intensity failure of the low-temperature exposure and the decreased 
absorption of the emulsion halide. Experimental methods of reducing 
the observed speed loss are just those to be expected on the basis of the 
above statement. ‘The wastage of an exposure by recombination appears - 
to take place to an equal extent at both room and low temperatures. 
Finally it is to be noted that the explanation given does not depend on 
the theory of high-intensity reciprocity failure ; in remains true whatever 

‘the explanation of the various features of a reciprocity failure curve. 
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XXIV. Note on the Angular Distribution of Radiation from Fairly Thin 
Targets Bombarded by High Energy Electrons 


By J. D. Lawson 
Atomic Energy Research Establishment, Harwell, Berks.* 


[Received December 19, 1951] 


At least four papers have appeared in recent years in which a calculation 
is made of the angular distribution of bremsstrahlung from targets 
sufficiently thin that ionization loss and straggling are negligible but — 
sufficiently thick that multiple scattering of the electrons is important 
(Schiff 1946, Lawson 1950, Kulenkampff, Scheer and Schittenhelm 
1951, Lanzl and Hanson 1951). In no two of these is the expression 
obtained for the angular distribution the same, nor have the method of 
derivation and the approximations used in the calculations always been 
clearly stated. This note is intended to clarify the position, by showing 
that there are two essentially different ways of carrying out the analysis, 
which lead to slightly different results. 

In the analysis which follows the small angle approximation will be used, 
and a gaussian law for multiple scattering will be assumed. The use of 
the small angle approximation means that the upper limit in integrations 
over all angles cannot be defined, for analytical convenience it will be 
taken as infinity. The following symbols will be used : 


EL =kinetic energy of electron. 

fe =rest energy of electron. 

t =thickness of target. 

f(@)  =angular distribution of the radiation from a target so thin 


that multiple scattering is negligible (this function has an 
angular width. of order y/H ; its analytical form will be given 
later). The function is normalized so that 


{p (0) 270d0=1. 


F(t, 0) =normalized angular distribution of scattered electrons after 
passing through a thickness ¢ of the target=(1/7ct) exp (—@/ct) 
where c is a parameter depending on the electron energy and 
target material. At high energies ¢ oc 1/H?. 

oot =fraction of electron energy lost by radiation in a distance 5¢. 

y =Luler’s constant=1-78. 


We shall now calculate the forward intensity and angular distribution by 
_two methods. The first is that employed by Schiff. The fraction of 


* Communicated by the Author. 
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energy radiated per unit solid angle in the forward direction from a target 
of thickness ¢ is clearly given by 


Ry=o[ fa F(t, 6) f(0) 2nOdedt, =. 2... (I) 
substituting for F(t, 0) and integrating with respect to ¢ 

Pee elec) | : —Bi(—€Jct) f(0) 2n6d8, . . . . (2) 
this may be simplified further if ¢ is not very small. Now f(@) is small 
when @>,./#, so that if 2/Hct is small then 62/ct will be small for values of 


@ for which f(6) is not small. The expression Hi(—é/ct) may then be 
replaced by the small angle approximation —In(y6/ct), and eqn. (2) becomes 


Se r 00 Hct pe 
Ry & [in (Stee) 10) 20 Bee te ie (3) 
=(a/ac) (In H’ctiu?—In dv). . . . . . (4a) 
where In A= | QO f(0) In (yH?6%fu2)d0.° 2 . 2... (46) 
0 


Substituting the approximate expression (Schiff 1951) 


iio) alae (PSMA) a a Se 3 (5) 
in eqn. (46) yields the remarkably simple result A= y, whence from 
eqn. (4 a) 

Pee (Gane anche), 2 ewer!) 2, mek Se 41 x ¢ (6) 


The assumption above that .7/H?ct is small clearly breaks down when ¢ is 
very small, nevertheless there is a useful range of ¢ over which eqn. (6) is a 
good approximation. . 

So far we have found only the forward radiation intensity. It is not 
possible to find a simple expression for the radiation at small angles,.but at 
large angles, where 0 is large compared with the angular width of (9), (i-e- 
p/H#) the angular distribution of radiation is simply given by the total 
distribution of electrons throughout the thickness ¢ of the target 


Reosyjp=(—e/n) Et (—@?/ct). . . ° . . (7) ; 

The value of R for small 6 may be found by judicious use of a French curve. 
The second method of calculating the angular distribution was first used 
in a crude form by Lawson and then refined by Lanzl and Hanson. The 


refined calculation will be given here. 
In this method we express f(#) as a sum of gaussian functions of different 


width, (in practice two components are probably all that can be justified), 
thus — 
f(0)= (a,/7) exp (—b,H?P}u?). . . . «.. (8) 
Now the angular distribution of radiation from a layer of thickness 6¢ at 
a depth t in the target is given by the convolution of (0) and F(t, @). Since 
Y2 


f 
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both these expressions consist of gaussian functions, which Conn ote 
simply by addition of their mean squared angles we get 


a,cot 0 
whence by proceeding to the limit and integrating from 0 to t 
o : i ; 6G 
tem 2 20,| 8 (—aoram) +8i(— Fam) |. - 00 


This expression is valid for all values of 0, and does not diverge for 0=0, 
and ¢ small. For ¢ and @ large it simplifies to eqn. (7). When 06=0 it 
simplifies to 

Ro=(e/mc) 2) a, in (1 =-cll*]b 7) ee eee ae 
In Lawson’s paper only one gaussian function was taken for f(#), and b, was 
taken as 2. 

It is interesting to note that if a single gaussian function is taken for f(6) 
and substituted in eqn. (4 6) we get A=b, the expressions (4 a) and (10 a) 
are then identical when ctH?/y?>1. 


Note added in proof.—Since this letter was written a paper by 
E. G. Muirhead, B. M. Spicer and H. Lichtblau has appeared (Proc. Phys. 
Soc. A, 65, 59, 1952), in which yet another formula for R, has been 


suggested. 
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XXV. The Specific Heats of Germanium and Grey Tin at Low Temperatures 


By R. W. Hitt and D. H. Parkrnson* 
Clarendon Laboratory, Oxford 


[Received November 16, 1951] 


ABSTRACT 


The specific heats of germanium and grey tin have been measured 
from liquid helium temperatures up to 170° and 110°x respectively. 
A comparison of the results with existing data shows that previous 
measurements on germanium do not represent the behaviour of the 
pure material, whereas those for grey tin are satisfactory. It is also 
shown that the temperature variations of the Debye characteristic 
temperatures of these elements are very large, and qualitatively similar 
to those of diamond and silicon. This behaviour is thought to be 
characteristic of lattices of diamond structure, but the possibility of 
the existence of anomalous contributions to the specific heat is also 
discussed. 

Values of the enthalpy, entropy and free energy of both substances 
have been calculated as a function of temperature and are presented in 
tabular form. 


§1. [INTRODUCTION 


For many years the Debye formula was regarded as giving a very 
satisfactory representation of the specific heats of regularly crystallizing 
solid bodies, although of course its theoretical shortcomings had been 
realized from the beginning. It is true that there were few substances 
where it fitted the data completely and it had therefore become the — 
practice to plot the Debye 6, as calculated from the individual measured 
points, as a function of temperature. Variations of the order of about 
10% were fairly common. 

However, the substances with diamond structure, particularly silicon 
and grey tin, showed very much bigger deviations and this was later 
shown to be the case for germanium also. It had therefore been proposed 
by Simon (1926) that one had to consider additional contributions to 
the specific heat in a way which will be discussed later in this paper. 

When it was later found by Blackman (1937) that the real frequency 
spectrum of a solid could differ very appreciably from Debye’s approxi- 
mation the question occurred whether the great deviations in the diamond 
type lattices could not have been explained simply in this way. The 
present work on germanium and grey tin was undertaken in order to 
SE a a capt tee ae Si A Scie ae Gait el ca ae i Se beret 
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extend the existing data down to lower temperatures. Also, in the case 
of germanium, which in the early experiments had been available only 
as a small and impure specimen, new measurements over a wide range 
of temperatures were thought desirable. 


§2. EXPERIMENTAL 


Two samples of germanium were used: the first was 35 g of Dupont 
‘ Hyperpure’ germanium in the form of granules of 2 to 3 mm diameter. 
This material is ‘ spectroscopically ’ pure, but may contain some nitrogen 
which does not show up in spectroscopic analysis. The second sample 
was a single ‘ button’ weighing about 10g, which was kindly supplied 
to us by Professor K. Lark-Horovitz of Purdue University. At the time 
of its arrival (September 1947) it was regarded as being of very high 
purity, having a specific resistance of 1-2 ohm cm compared with the 
theoretical value of 60 ohm cm for the ideally pure material. 

For the work on grey tin two samples were also available. The first 
was 259 of rather coarse powder believed to have a total impurity 
content of about 0:05°%. The second sample was kindly supplied by the 
Tin Research Institute, who believed it to be of the purest tin available. 
It is their opinion that existing methods of estimating impurities in tin 
are unsatisfactory, and that no significant analysis can be given. The 
sample used was 100 g, and was similar in form to the first sample. 

The measurements on both samples of germanium and the first sample 
of grey tin were made in an apparatus which has been described elsewhere 
(Parkinson, Simon and Spedding 1951). In this apparatus, the lowest 
temperatures were attained by pumping away liquid helium which had 
been condensed into the calorimeter. This procedure was found to be 
unsatisfactory for grey tin, since the helium was quite strongly adsorbed 
on the powder sample, giving errors of up to 20° due to the progressive 
desorption of helium as the sample warmed up. A somewhat different 
calorimeter was therefore used for the low temperature measurements 
on grey tin, differing from the other in that a separate chamber was 
provided for the liquid helium. It was found necessary to admit about 
10° moles of helium exchange gas into the calorimeter in order to establish 
thermal equilibrium between the sample and the thermometer. A com- 
parison with the amounts adsorbed in the other calorimeter, together 
with the approximate error introduced by it shows the error due to 
adsorption in the second calorimeter to be less than 1° above 7° K. 

The temperature ranges covered in this work were as follows : germanium 
(both samples) 4 to 175° kK ; grey tin (first sample) 12 to 110° kK; grey tin 
(second sample) 2 to 20° K. 

Above about 10° x, the results obtained are believed to be accurate to 
about 1%. At lower temperatures, the accuracy of measurement decreases, 
due to the fact that the heat capacity of the specimens becomes very 
small and comparable with that of the empty calorimeter. Under these 
circumstances, the unfavourable ratio of surface area to volume of the 
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sample, which is inherent in work with granular or powder specimens, 
adversely affects the thermal insulation which can be achieved, and with 
it the accuracy of the result. 

In the case of grey tin, a further limiting factor is the adsorption of 
exchange gas by the specimen at helium temperatures. The effect of 
this is difficult to estimate, but would appear to be rather less important 
than the effect of poor thermal insulation, yet the final accuracy may be 
as low as 20% below about 4°x. Thus these measurements do no more 
than establish the absence of marked anomalous behaviour between 
2 and 4° x. 


Fig. 1 


20 40 60 80 100 120 140 {60 


Oo 5 Te) 15 20 25 + 30 35 


The specific heat of germanium. The upper curve shows the results at higher 
temperatures. 
Abscissa: temperature ° K. Ordinate : C> cal/deg. g atom. 


§3. RESULTS 
(i) Germanium 

The results obtained with both samples of germanium are shown in 
fig. 1, and it may be seen that both sets of results lie on the same smooth 
curve. These results are very different from those obtained by Cristescu 
and Simon (1934); no peak is found in the curve and the specific heat 
values are everywhere smaller. Over a wide range of temperature this 
difference is as large as 50%. 

Smoothed values of the present results are given in table 1 together 
with values of the specific heat at constant volume and the Debye 0 
deduced from them. Also included in this table are values for the entropy 
iS, the enthalpy H and the Gibb’s free energy G. The correction from — 
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C,, to C,, could only be made approximately since the expansion coefficient 
and compressibility are not known. Use was therefore made of the 
formula due to Nernst and Lindemann 


C,—C,=0-02140,77/T,, 
where 7’, is the melting point of the substance. This formula is, of 
course, only approximate, but the corrections applied nowhere exceed. 
1-5% and are negligible below 70° kK 


Table 1. The Specific Heats and Thermodynamic Functions of Germanium 


2-59 
2-92 
3°22 
3°77 
4-26 
4-68 


(The thermal unit is the calorie per g atom throughout.) 


(ii) Grey Tin 

The data of Lange (1924) extends down to a little below 20°x 
measurements were therefore made in the temperature range 2 to 30° K 
Results obtained with both specimens agree very well with Lange’s values 
in the region where these overlap. Further determinations in the 
temperature range 60 to 100° x confirm this agreement. The present 
results are shown in fig. 2 together with those of Lange; both sets of 
results have been used in compiling table 2 which contains smoothed 
values of C,,, C,, 0, S, H and G. Corrections from C, to C,, were made 
in the manner described above, the difference being 2°% at 100° xK and 
negligible below 40° kK 


§4, Discussion or REsuuts 


The results described above show that the values of, 6 for germanium 
and grey tin vary with temperature in the same way as those for diamond 
and silicon. In all cases, 0 tends to a constant value 6,, at the higher 
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temperatures, passing through a minimum at lower temperatures. If 
6/8 ., is plotted against T'/@., as in fig. 3, this essential similarity is clearly 
displayed. It is seen that the minimum becomes deeper and occurs at 


Fig. 2 


~ ao 10 15 20 
The specific heat of grey tin. The upper curve shows the results at higher 
temperatures. 
<x =Lange,- ©=this research. 
Abscissa : temperature ° K. Ordinate: C, cal/deg. g atom. 


Table 2. The Specific Heats and Thermodynamic Functions of Grey Tin 


(The thermal unit is the calorie per g atom throughout.) 
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lower temperatures as the atomic weight increases. The data for silicon 
used in this diagram is that of Keesom and Pearlman (1951 a); this 
data, together with some for germanium, has recently been summarized 
(Keesom and Pearlman 1951b). It would appear that this data for 
germanium is in good agreement with that presented here. 

In the case of diamond, the minimum has been explained purely in terms 
of lattice theory by Smith (1948) and a recent note by Hsieh (1951) 
reports that a similar calculation for germanium yields values which are 
in reasonable agreement with the present results. It thus appears probable 
that a similar treatment would account for the behaviour of silicon and 
grey tin also. 


Fig. 3 


O on 0-2 O-3 
The reduced Debye temperature as a function of reduced temperature for the 


diamond structure elements. 
Abscissa : 7/0,,. Ordinate: 0/0,. 


An alternative approach to the problem is that proposed by Simon 
(1926), in which it is postulated that the deviations from the Debye 
theory are due to an additional contribution to the specific heat arising 
from some excitation process. Such a process could contribute a term 
to the specific heat which would be represented by a Schottky function. 
This view gains support from the fact that the experimental curve can 
be fitted quite accurately by the sum of a Debye and a Schottky function, 
using the characteristic temperatures given in table 3. Perhaps too much 
importance should not be attached to this fact, since the effective 
vibrational frequency spectrum may be rather similar in both cases, and 
the specific heat is insensitive to the form of this spectrum except at 


, of Germanium and Grey Tin at Low Temperatures 315 


the lowest temperatures (Katz 1951). It must also be remembered that 
there is no known origin for a suitable excitation process in the pure 
materials at these temperatures ; indeed there is no need to invoke such 
& process at all for diamond. On the other hand, the existence. of a 
negative expansion coefficient for silicon (Simon and Bergmann 1930, 
Erfling 1942) would seem to indicate that the behaviour of these substances 
is not perfectly understood. The expansion coefficients of silicon and 
germanium are to be measured shortly in this laboratory, and the effects 
of impurities investigated. 

A further matter requiring explanation is the very large discrepancy 
between the present results on germanium and those of Cristescu and 
Simon. The specimen -used by the older authors was small, but the 
difference between the two sets of results seems far too large to be 
explained in terms of experimental error. Aigrain and Dugas (1950) 
have shown that the heat content of the Cristescu anomaly can be 


Table 3. The Debye and Schottky Characteristic Temperatures for the 
Diamond Structure Elements 


Element 


Diamond 


Silicon 
Germanium 
Grey tin 


explained in terms of a concentration of impurities of the order of 
1019 impurity atoms/cm*. If the remaining difference between the two 
sets of results is attributed to impurities in the specimen of the older 
authors, the specimen must have been very impure indeed. While this 
may have been tlie case, the entropy difference at 200° K is of the order 
of RIn2; it would therefore appear that about one extra degree of 
freedom per atom is required in order to explain the discrepancy, and it 
is difficult to see how this could be attributed to impurities. A final 
- decision on this point must await experiments carried out with controlled 
amounts of impurity ; such experiments would be of very great interest 
in the study of all the diamond structure elements. 
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XXVI. The Relation between Thermo-Optic and Piezo-Optic 
Phenomena in Crystals 
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SUMMARY 


It is pointed out that the change in refractive index with temperature 
of a crystal is different from what is calculated from the accompanying 
change in volume and the piezo-optic coefficients. ‘The difference, which 
is a pure temperature effect, is explained as being due to the change in 
polarizability of the atoms produced by a change in the amplitude of 
vibration. The polarizability («) can be expanded as a Taylor series in 
the changes of the distance (r) between the atoms and it is found that 
while the piezo-optic coefficient depends only on d«/dr, the pure temperature 
effect is a function of 0?«/dr?._ Making use of the experimental data, 
the values of « and its first two derivatives can be determined. These 
values are found to be of the same order as those deduced from the 
intensities of Rayleigh and Raman scattering of light. The theory 
predicts that dn/dT should vary as the coefficient of cubical expansion 
at different temperatures and this is verified to be true. Finally, 
calculations are made of the thermo- and piezo-optic coefficients, 
considering the electrostatic interaction between the atoms. These do 
not adequately explain the observed facts, since no provision is made for 
the distortion of electron atmospheres around the atoms and_ the 
consequent changes in polarizability. 


§1. INTRODUCTION 


WHEN a transparent solid is subjected to a change of temperature or 
pressure, its refractive index, n is found to change. The two quantities 
dn/dT' + and dn/dp may be referred to respectively as the thermo-optic 
and piezo-optic coefficients. Measurements of the former exist for a 
number of crystals over a range of wave-lengths and temperatures (e.g. 
Micheli 1902, Reed 1898, Ramachandran 1947, Radhakrishnan 1950, 
1951). The piezo-optic coefficient has not been measured directly for 
any crystal, but it can be calculated from the piezo-optic constants ¢,;, 
or the elasto-optic constants p;; of the crystal. These latter have been 
measured for many crystals (e.g. Pockels 1902, Ramachandran 1950, 


Burstein and Smith 1948, 1949). The present paper is a study of the 
See Catron ee ee ee 
* Communicated by Professor R. 8S. Krishnan. pig 
+ As against ¢ for temperature, used by the authors in earlier. publications, 


T is used here in accord with thermodynamic symbolism. 
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relation between dn/d7 and dn/dp. Both on the application of a 
compressive stress and on reducing the temperature, the specific volume 
of a substance diminishes. If the changes in volume are made the same 
in both the cases, are the changes in refractive index also the same ? 
Calculations show that they are not; in fact, in some crystals, even the 
signs of the refractive index change are opposite in the two cases. It is 
therefore clear that a factor affecting the refractive index is operative, 
purely caused by a change in temperature, in addition to that due to 
thermal dilatation. The existence of this factor has been discovered long 
ago (Pockels 1902); but no explanation appears to have been given for 
it. An explanation is attempted in this paper, in terms of the polariza- 
bilities of atoms and the influence of thermal vibrations on polarizability. 


§2. CoMPARISON OF dn/dT’ anpd dn/dp 
The dispersion formula for a crystal can be written in the form 


t(n)=40 2 aes 


where J, is the number per unit-volume of the r type of oscillator which 
has a polarizability «,. The form of the function f(n) depends on the 
so-called ‘ internal field ’ in the crystal. Two forms are commonly used, 
viz., the Sellmeier form f(n)=n?—1 and the Lorentz—Lorenz form 
J (n)=(n?—1)/(n?+-2). In a paper appearing elsewhere, the authors have 
shown from the experimental refraction and absorption data that the 
first form is more applicable to crystals, which is therefore adopted in 
this paper. Replacing the summations on the right-hand side of eqn. (1) 
by a term Na, where N is the total number of oscillators per cm? and 
« the average polarizability of each, we have 


n?—1=47rNa. on bine <i eaek a ae an 


For an analysis of thermo-optic and piezo-optic data in the visible region, 
where the absorption frequencies are far away, we can justifiably use this 
simplified ‘ one-term dispersion formula ’. 

When p or 7 is changed, both N and « will alter. The change in N 
is completely determined by the change in the density, i.e., 


aN _aN dp dN _ aN dp 

dp dp dp’ aT dp aT’ . 5 . . Py (3) 
and since N is proportional to p, dn/dp=N/p. The polarizability « is, 
in general, a function of both the density and the temperature i.e., 


da A. dp du 0x dp (du 
dp dp' dp’ af ap'at (an). (4) 
Now 
dp dp 
Aron and ipa Pre Te Na ope pee (5) 


where 6 and y are respectively the compressibility and the eee 
of cubical expansion. 
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Differentiating eqn. (2) and substituting from eqns. (3), (4) and (5) 


we have 
dn ,n?—1 p (dx F 
ee ela ee 


dn n2—] p (Ox “1 (0a 
a= 7a [+2 (55), +a(ar),| 


It is thus possible to separate dn/dT into three parts, due to (a) the 
change in number of oscillators per em3, (b) the change in polarizability 
caused by thermal expansion, (c) a pure temperature effect on polariza- 
bility. These parts are respectively denoted by P, Q, R in conformity 
with the earlier notation of one of the authors (Ramachandran 1947). 
We have from eqns. (6) and (7) 


Pp on aN n2—1 
~\ON). aT "7 on 


j= On\ Oa dp _ ldn n?—1 
Ox), Cf.ar | 375 = |: Gren A) 


pa (er (d\ _ dm. yan 
= (G) (an) ,= an Bap 
These equations enable us to calculate P, Q, R in terms of the macro- 
scopically measurable quantities dn/dp, dn/dT, y and B. It is also evident 
that 

d B 


n dn ; 
apuPt+OteR; pain ams ae cs 1 (9) 


Using the data for y and £ available in the standard tables, values of 
P,Q, R have been computed for a number of crystals and are shown in 
table 1. 


and 


? 


Table 1. Values of P, Q, R in Crystals at 5893 A 


e 
Crystal nm 10°dn/dT 105, (=) LOO fle (Qin | 


Op T 

NaCl 1544 =3-8 12-0 (i 24uee =) ame 20-08) 2-8 
KCl 1490 —36 11-4 0-23. —4-7 Dl Pe 1-Oaeo | 
KBr 1559 .—3-6 12-0 0:35 —5-6 Td . 10-6 1 24 
KI 1-667 —5-0 13-5 0-430 = 7-2 4e, #0818 
LiF 1-392)¢ =1-3° ° 10-2 Oren erm 2-4 28.8 8H 
CaF, LAsieeelD 5-7 0-259) 22:1 O27 e202 e385 
MgO 1-736. 11-6 8.9 1 0.4) 1s 8 ey GE 
Diamond 2-417 +41-0 1:3 —0:28 —1-3 166 -+0:64 2-7 
Quartz w 1544 —0-65 35 (Seer Game 0-8. 0-6 205 

feel 553 8 e=0-Th) 3 C40Me Ga 0:20 220-7 03 
Topaz « 1596 +0-81 1-52 0-03 —06 0:55 +09 06 

6 1610 +0:90 1-52 O04 0-7. 0-60, 0-99 80-7 

y 1:6778 +0:80 1:52 0-04 —0-7 (-Geyaeeoo. Os 
Calcite. wid-658 .---0-07 158 075 —0-8-° —0-06 41-37 00-05 
| 21486 (1:06 158 —013° —06 045 =+0-9. 05 
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The table shows that, irrespective of the sign of dn/d7’, Q is always 
positive and R is generally. positive. Thus the polarizability mereases 
with expansion and usually increases with temperature. The universal 
negative value of P is clearly due to the decrease in density on heating. 

The table has been so arranged that the first half consists of crystals 
which are usually classed as ionic and the second half consists of crystals _ 
having covalent linkages. Calcite is considered covalent since its 
optical properties are determined by the CO, ion, which in itself is covalent 
in character. It is seen that in ionic crystals the magnitude of the pure 
temperature coefficient R is smaller relatively to Q@, than in covalent 
crystals. This is only a general trend, diamond being a notable exception. 
Further our derivation of P, Q, R pertains, strictly speaking, only to 
isotropic crystals, so that the values for calcite should not be regarded 
as accurate. Finally, it should be mentioned that the experimental 
determination of elasto-optic constants are not very accurate. Since 
the magnitudes of P, Q and RF are very much affected by the value of 
dn|dp used, the whole of table 1 should be regarded only as semiquantitative 
in character. 


§3. EXPLANATION OF THE PURE TEMPERATURE EFFECT 


Consider the simple example of a diatomic molecule, and let 7, be 
the equilibrium inter-atomic distance at the temperature 7’. As the 
temperature is raised, 7) increases and further the instantaneous separation 
distance 7 executes vibrations of increased amplitude about the mean - 
value ry. If v, be the frequency of vibration and A its amplitude, then 
the displacement of an atom at any instant from its equilibrium position 
with respect to the other atom may be written 


Ar=A cos 271,t. oi) 20> ee eae Se a 


Now the optical polarizability « may be expanded in a Taylor series 
in Ar as 


0 1 /e? 
a(7% + Ar)=a(7%9)-+ @) Ar+ = (53) (Ary?+..... s coe ia 
“ 0 


We shall hereafter neglect higher terms than the second since Ar<r. 
We have assumed that the molecule vibrates harmonically about an 
equilibrium position, since the presence of non-linear restoring forces 
would only result in the mean distance being different from 7) and in the 
presence of higher harmonics of v,. These will not affect our discussion. 

Substituting value of dr from eqn. (10) in eqn. (11) and integrating 
over a cycle of v,, one finds the mean polarizability to be 


A? (a 
a=alrolt | (Fa) - risk dace Sen) 


This picture is valid for a diatomic molecule. But as a simplified 
extrapolation, we suppose it to be valid for a crystal also. In doing 
this we neglect many finer details. For instance a crystal will have not 


between Thermo-Optic and Piezo-Optic Phenomena in Crystals 321 


one, but many thermal lattice vibration frequencies. We shall provide 


in a rough way for this by taking A to be the root mean square amplitude 
of the thermal vibrations. From eqn. (12) we have 


(5 _ [0a dry 02x 0A?/4 ; 
aie (aban + a), a Jeo = 08) 


When the crystal is compressed isothermally, all linear distances are 
altered, and since 7, 0c p~ 18, 


Gin G2 1G] Opens, Lae tare re ne a (14) 


Further, the vibration amplitude A is, at least to a first order, only a 
function of the temperature of the crystal, so that we may write 


(0A?/dp)p=0. order, | yt een wet LD) 
Hence we obtain 
Ox al Moles ; 
Likewise 
(elca 1 dA? (07a 
(jn).= iar (aa),: . . ° . . . (17) 


Thus we see that the pure temperature effect arises from the presence 
of the additional factor dA?/d7’ governing thermo-optic phenomena. 
Since this factor is always positive, the sign of R depends on (0?«/dr?), 
i.e., the coefficient of the quadratic term in the Taylor expansion. 

From the above equations it is possible to find the relative orders of 
magnitude of «, (dx/0r), and (0?a/dr?)). We have 


Ue iP 
. . . . . . . 8 


imme coe] (Cie 8 aA t/dt © (19) 
(=), ae oe 


From table 1, it will be seen that P/Q and Q/# are of the order of unity. 
Now the following values may be taken as representative of the various 
quantities involved: r=2x10-§ cm, y=104, A=0-2 a (Lonsdale 1 948) 
. at room temperature. Neglecting zero-point energy, we may estimate 
dA2/dT simply as A?/T7', T being 300° x. From these values one obtains 


a ~ 10-8(8a/Ar) ~ 10-16(22x/Ar2). 2 2 2 2. (20) 


and 


Thus, if the distance 7 is measured in Angstroms, « and its derivatives 
are all of the same order of magnitude. Substituting this result in the 
Taylor expansion (11), and remembering that Ar is of the order of 0-2 a, 
we see that terms in the Taylor expansion beyond the second may be 
safely neglected. 


SER. 7, VOL. 43, NO. 338.—MARCH 1952 Z 
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$4. ComPaRISON wiTH Licut ScaTTERING Data 


It is well known that the intensities of Rayleigh and Raman scattering 
by a transparent medium are also governed by the quantities «, da/dr 
and @a/dr2. A rigorous theoretical treatment which enables one to 
predict the scattered intensities in any particular crystal is a very 
complicated problem, and is in any case beyond the scope of this paper. 
But as a very rough approximation one may take the well-known 
expressions for the intensity scattered by molecules (cf. Bhagavantham’s 
‘Scattering of Light and the Raman Effect ’, Ch. V) and modify these, 
to suit the case of crystals. Using the Sellmeier dispersion formula, one 
obtains the following expressions for the ratios between Rayleigh scattering 
intensity Ip, first order Raman scattering J, and second order Raman 
scattering I, : 


In, (P+Q\*NETB : 
ifr (5 ys | oO (sam) (21) 
and 
I, 4 da/or \2 99 
i= 3(seg) - ao” SP See 


where N is the number of oscillators per cm’, and o is the amplitude of 
oscillation, given by 


Dart 2 RY cs aio: ee ke a an cee 


yw being the reduced mass of the oscillator. 

These expressions are not rigorous formulae for scattered intensities, 
but demonstrate how the ratio of « to its derivatives influences the 
intensity in the three types of scattering. Experimental data on 
scattered intensities are very scanty. But from values in diamond 
(Krishnan and Narayanan 1950) and quartz (Krishnan and Chandrasek- 
haran 1950) and from the experience of workers in this laboratory 
(Chandrasekharan, personal communication) one concludes that, in 
crystals, first order Stokes-Raman lines have an intensity lying between 
0-1 to 10 times the intensity of unmodified scattering, while the second 
order lines have a, few hundredths of the intensity of the first order. 
Taking 8 as 3x 10~" and N as 10?5, which are representative values for 
many crystals, and (P+Q)/P as 2, we can evaluate the ratio of « to its 
derivatives from eqns. (21) and (22). We then get the results 


a/(Ax/Or) ~ 108 and (dx/dr)/(d2x/Ar2) ~ 108,  % 4(24) 


which are in accord with the orders of magnitude obtained from’ thermo- 
optic and piezo-optic data. 


§5. THe VartaTion or dn/dT wrrn TEMPERATURE 
Substituting eqns. (16) and (17) in eqn. (7) we have 


w—l1r, (du v—1dA2 1 (a 
¢ 23 —_—_—_—"—_——- alll < == Ms 2A 
ea 2n (5), : R= Qn aT 4a. = (Se), eee 
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Since dA?/dT is proportional to the rate of change of thermal energy with 
temperature, i.e., O,, we may write 


dn 
Gp re R=ay 0, 2 2. es (26) 


where a and 6 are practically constants. Now for solids, C,=C, and 
further, C, is approximately proportional to y (Gruneisen’s law). Hence 
we get the approximate relation ; 


dn 
ap © Cp ORG ROC py tes vet Dhty Be Peee Cour) 


This in fact is found to be the case in crystals, dn/dT’ being very small 
at low temperatures and increasing with temperature in a manner similar 
to y. Of course the proportionality cannot be expected to be strictly 
valid at all temperatures, since we have no reason to assume that da/dr 
or 0?x/dr? is unaffected by temperature changes. But in diamond and 
quartz it holds to a fair degree of accuracy as shown in table 2. 


Table 2. Relation between dn/dT and y in Diamond and Quartz 


Diamond Quartz 
Temp. 
(deg: c) . (1) (2) (3) (4) 
4105 dn/dT 10°, (dnjdT)/y | —10°dnjdT 105,  (dn/d7)/y 

—150 . 0-01 — — 4-4 1-9 2°3 
—100 0:19 0-09 2-1 4:8 2-52 1-9 
—50 0:50 . 0-17 2-9 5:2 2-98 1:8 
0) 0:84 0:29 2-9 5:65 3°36 1:7 
50 1-19 0:43 ° 2°8 6°15 3°68 1:7 
100 1-54 0:55 2°8 6°75 4:00 1:7 
150 1-91 0-65 2-9 7-40 4:33 1-7 
200 * 2-41 0:77 3-1 8-20 4-66 1-8 
250 2-67 0-85 3-1 9-20 5:04 18 
300 2-98 0:92 3:2 10-5 5:49 1:9 
350 © 3°32 0:98 3-4 13-0 — — 
400 3°62 1:03 3:5 14-6 6°74 2-2 


(1) Ramachandran, 1947; (2) Krishnan, 1946; (3) Radhakrishnan, 1947 ; 
(4) Sosman, 1927. 


For a crystal having a dispersion formula of the form 

Fea hi (veer ee Be ke (28) 
the variation of dn[dT with wavelength can be expressed by the formula 
(Ramachandran 1947) 
) dn n*?—1 byo* BO 
dP ton neve ‘ae 

where y,7 is defined by the relation 

L dry 1 dy 
Le eT EME 3 GWU 
Z2 


(30) 
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It is found that yp is in general positive, is very small at low temperatures, 
and increases with temperature. Since both the terms on the left hand 
side are proportional to y, one should also expect that 


Xr © ¥. (31) 


This enables one to understand the trend of the observed variation of 
x7 with temperature. 


§6. CancuLaTion oF dn/d7T’ anp dnj/dp vor a Cusic Latrice 


In crystals like the alkali halides, the polarizability and hence the 
refractive index is mostly due to one type of ion. Let «, be its polariza- 
bility in isolation (e.g. in the gaseous state) and «, its polarizability in 
the crystal, which has only one component for a cubic crystal. Then 
x, is given by (Bragg 1924) 

3 cos? 6;—1 


re ae ia En 


p= OH gf 2! g? 

J 
where 7; is the length of the vector r; joining an atom j (coordinates 
2X, Y;, 2;) to a chosen origin and @, is the angle between r; and the direction 
of the incident electric vector. But the last term in eqn. (32) is identically 
equal to zero for a cubic arrangement of ions about a central ion* (Mott 
and Gurney 1940). Thus 


eT ere Shoe AG A ce RS 


so that 
Ou 0 
p (),= Pee Tre SN RES 
and from eqn. (6), 
dn ,n?—] 
at oI Rewer re ee eT 


The right hand side is proportional to P and the theory leads to the 
result that Q@=0. This is not in agreement with experiment, since Q is 
finite and positive for most crystals. This is not surprising since the 
theory takes into account only the purely electrostatic interactions 
between neighbouring ions. It is conceivable that the polarizability 
would alter owing to deformation of the electron atmosphere of the ions 
when the crystal is compressed and the effect of this has been neglected. 
Actually, even the result (33) is contrary to facts for the values of «, 
and «,, derived by Burstein and Smith (1948) from the experimental data, 
are not equal to each other. The difference arises because one cannot 
assume a priori that the ion has exactly the same structure in the gaseous" 
and the crystalline states. Owing to the packing in the crystal, the 
electron atmosphere of the ion would suffer a certain amount of distortion, 
leading to a change in its inherent polarizability. 


* We are grateful to the referee for certain suggestions on this point. 
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On the basis of the same considerations, one can also work out the 
value of the pure temperature effect, R. As a simplification we may 
picture the random thermal motion of the atoms in a cubic vend to be 
as follows : A particular atom is considered to be executing simple 
‘harmonic vibrations with all possible orientations, the other atoms 
being considered to be at rest. If Ar is the amplitude of vibration for a 
particular orientation, with components Aa, Ay, Az, then its time-averaged 
contribution to the polarizability is, to a second order, given by the 
generalization of eqn. (11) as 


meet | (ara 02x . o2 
sent) srs) (8) 0 


Writing 
Ax=Arsingdcos#, Ay=Arsindsind, Az=Arcosd, (37) 


and integrating over all possible values of 6 and ¢, one obtains the mean 
value of (4a)? over all orientations as 


: [ [7 (4a)? 2 sin dg do — ee 35) 


4 }o Jo 3 


and the same values for (Jy)? and (4z)?.. The mean values of the cross- 
products Ax Ay etc., vanish. Thus 


(Ar)? (a Ba  ay\ (Ar)? 


Aa= 


Now, 
V2 {(3 cos? 6,—1)/rF}=0, . . . . «. (40) 


for every atomj. Thus, even if the summation in the last term in eqn. (32) 
does not lead to zero, 4x,—0, provided the vibrations of the individual 
atoms are isotropic. In any case, we obtain the result that (d«/0T'),=0. 
On the other hand, the fact that the coefficient FR is finite in cubic crystals 
shows that actually (d«/07'), is not equal to zero. As in the case of the 
' piezo-optic coefficient, this discrepancy between. theory and experiment 
is to be attributed to the fact that we have neglected the distortion of 
the electron atmospheres caused by thermal vibration. Such effects 
appear to be quite sensible, judged from the magnitude of the discrepancy 
both for Q@ and R. It follows therefore that one cannot expect to furnish 
a reasonable theory of photo-elasticity in crystals, based only on con- 
siderations of electrostatic interaction between neighbours. 


‘We are grateful to Prof. R. S. Krishnan for his interest in this work. 
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SUMMARY. 


A method is described for calculating the energy of the lowest state 
% of a divalent metal. This approximates to the well-known cellular 
method, but assumptions made regarding y,, together with the use of a 
single empirical ionization potential, render unnecessary the explicit 
calculation of an ion-core field. An expression for the total energy of 
the valence electrons is obtained, the free-electron approximation being 
used tentatively in evaluating all terms except the energy of the lowest 
state, and this is then applied to the calculation of the atomic radii, 
compressibilities, and cohesive energies of the divalent metals. The 
change of volume for a range of pressures up to 100 000 kg/cm? is also 
calculated and found to be in good Sareea with the experimental 
data of Bridgman. 


$1. INTRODUCTION 


Tue difficulties encountered when an ion-core field of the Hartree type 
' is used in calculating the cohesive energy of a polyvalent metal have 
been demonstrated in a recent application of the cellular method to 
magnesium (Raimes 1950). In particular, it is necessary to know the 
correlation energy of the valence electrons in a free atom, which, at 
_ present, can only be estimated in a heuristic manner. ‘This difficulty 
might be overcome by the use of a semi-empirical ion-core field, like that 
of Prokofjew (1929) for sodium, but so far no such field has been constructed. 
for a polyvalent metal. The present work attempts to circumvent both 
the drawbacks of the Hartree field and the labour of constructing a semi- 
empirical field by extending to polyvalent metals an approximate method 
first applied by Frohlich (1937) to the alkali metals. Only the divalent 
metals are considered here, because it is believed that the various approxi- 
mations involved apply best to these metals, but the treatment can easily 
be made quite general. 

Although the method was originally designed for the purpose stated 
above, its greatest advantage has been found in calculating changes of 
volume with pressure. The total energy of a metal is here obtained as 
an analytical function of the atomic radius, so that the pressure corres- 
ponding to a given volume can be computed directly with an accuracy 


* Communicated by Professor H. Jones. 
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greater than in the ordinary cellular method, where the total energy is 
obtained as a tabulated function. It will be seen in §4 that the 
theoretical curves giving compression as a function of pressure, in the 
calculation of which no empirical data for the solids have been used, 
are in good agreement with the experimental results of Bridgman (1941, 
1948), for a range of pressures up to 100 000 kg/cm’. 

The possibility of employing the present method depends upon a result 
of the work on magnesium mentioned previously. The field there used 
in computing yy, the wave function for the state of lowest energy, consisted 
of two parts—the field of the ion-core and the field of the valence electrons. 
Two alternatives were considered for the latter part, namely, the field 
due to all the valence electrons and that due to electrons with one type of 
spin only, and it was found that the calculated total energy of the metal 
was not significantly different in the two cases. This suggests that pp 
could, with slight loss in accuracy, be calculated using the ion-core field 
only, as for the monovalent metals, provided that the Coulomb interaction 
of the valence electrons is accounted for by an additional term in the 
expression for the total energy. The Coulomb term, per atom of the 
metal, is just the self-potential energy of the valence-electron charge 
distribution within an atomic sphere. 


§2. THe State oF Lowrst ENERGY- 


According, to the argument presented above, a good approximation 
to the wave function (7) of the lowest state may be obtained by solving 
the Schrédinger equation 


dh, 2 dehy | 
ate eV Oe 


subject to the boundary condition 


(7), =° ‘Se 


where V (7) is the potential energy of an electron in the field of the ion-core, 
and 7, is the atomic radius.” 


We wish to evaluate ¢9(r,) for atomic radii in the neighbourhood of : 


1%, where ry is the value of r, for which ¢, reaches its minimum. It will 
be assumed that, in this range of 7,, %(r) is‘nearly constant over most 
of the atomic volume. Furthermore, we shall confine ourselves to the 
region lying outside the ion-core, so that it is sufficiently accurate to put 
V(r)=—4/r, for a divalent metal ; it is, of course, necessary that 7, lie in 
this region, and it will be seen later (table 1) that this is so for all the 
divalent metals except mercury. The calculation may follow almost 


* In this paper, unless otherwise stated, energies are in Rydberg units and 
lengths in Bohr units. 
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exactly that of Fréblich (1937) for the alkali metals, but the following 
modification has the advantage of showing more clearly what approxima- 
tions are essential to each stage of the work. 
Writing R=rj,, eqn. (1) becomes 
a?k 
Tee + [e9(r,)— V(r) |R=0, ak a eee Meese (3) 


and the boundary condition (2) becomes 
dR R 
. ae Ee Gea aie 
It is to be noted that, through eqn. (4), R is a function of r, as well as of r. 
Hence, denoting differentiation with respect to r, by a dash, we have 


(4) 


d?R’ Eee 
qe leo V(r) Me’ + <9’ R=0, The pay ee 
or, using eqn. (3), 
UL MTs : 
Therefore, 
ts ts d?R dese adn -_ dik; 
, 2 “ta i ekg, Sy SEDER T SS. = (ie ek ts Doe ie 38 
a [ Redr=["| RTS Roe |ar=[2B'S, oO 
We now differentiate boundary condition (4) with respect to r,, and 
obtain 
PR dk’ Lae 
dr is dr yor, = alee ee | 
so that eqn. (7) reduces to 
mics d?R 
cof Redr=| RFE | =(M el ROMP. (9) 


So far the treatment (cf. Bardeen 1938) has been exact, but now, 
consistent with our assumption that %, is nearly constant, we make the 
approximation that 

4a 4a : 

= rele P= % TARO P=1, 
when y, is normalized to unity in the atomic sphere; in other words, 
we assume that the charge density at the surface of the atomic sphere 
is equal to the mean charge density taken over the volume of the sphere. 
Eqn. (9) now becomes 


Te CU (gy eeen eee i cece nae (LO) 


Using the further approximation that V(r,)=—4/r, in the required range 
of r,, eqn. (10) may be integrated to give. 
, cs au ree, Realy 


€9(7)= rs aa r, F) 


where A is a constant. 
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Since ry is the atomic radius corresponding to the minimum eigenvalue, 
it follows that A=2r,? and 


cote)=2 (35 — =). Se 8 85, 


This differs from the expression obtained for the alkali metals only in 
the factor 2. It remains to determine 7. 
Let po, be the wave function for atomic radius 79, so that 


Ih,9 , 240, 0 


ate tea colo) — VF) Mo, 00) es 
and | 
dbo, o 
ae =0. Pe ei ae nll tes tne 
( dr yea : ( = 
Also, let us write 7) and its corresponding eigenvalue ¢,(7,) as 
boil. af (Cue 9 aes ns 
and Ege?) ==<ol%o) 12705) ee ee 


Substituting (15) and (16) in eqn. (1), and using eqn. (13), we find 


Pf 2 dboodf , 24f d 
i, dry deirdre.) eee 


We shall only be interested in the values of ds)/dr and yy at r=ry, where 
dip, o/dr is zero, so that, for our purposes, eqn. (17) reduces almost exactly 
to 


fee , 
oe ae nf=0. 2 eo) ee eee 
Hence f()=B — Bt +C ee : 


where f= ,/y, and B and C are functions of 7,. Now, when 8 > 0, f(r) 
must tend to the constant value unity, so that, provided f is small, it is 
reasonable to neglect the cosine term and write, near r=/r, 

sin Br 


ie Bras ev) gee (19) 
From eqns. (14), (15), and (19) we obtain 


1 diby 
(5; Ake 3 ¥ (Bro cot Bro—1), a © ee (20) 


which will be required in the evaluation of 7). It may be observed that 
in deriving this expression it has not been necessary to determine bo, 0 
explicitly, as in Fréhlich’s treatment, nor has any assumption been made 
regarding V(r); in fact, the only real approximation has been the neglect 
of the cosine term in f(r) 

When r, tends to infinity, eqn. (1) becomes that of the single remaining 
valence electron in a free atom when one valence electron has been 
removed, and —e,(7,) is then just the second ionization potential of the 
neutral atom, which we shall denote by J,. In the case of the alkali 
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metals, of course, the first ionization potential is used, but the rest of 
the procedure follows closely that of Fréhlich. In the region where 
V(r) is hydrogenic the Schrédinger equation for the state 7% of the valence 
electron in the free, singly-charged, ion is 

dys 2dr 4 

Ft ae (; 1.) y=0. Me re ee OL 
This equation may be solved to give an asymptotic series for ws, of which 
only the first two terms need be retained ; thus, 


pase ryt E | ; 3 Pore : (22) 


2wr 


where w?=TJ,, and m=2/w. 
Now ¢ is simply that particular function %s, which corresponds to the 
value of 8 given by 
+ 


Bn2= —€q(19) —La= re Be his an. tema ae ML 8.) 
0 
Hence, from eqn. (20), we may write 
ld 1 
€ i) ia; (Bo) Cot Boto—1), 


and, using eqns. (22) and (23), the resulting equation for 7, is 
2+w (2—w)? %  2—-—@ 

2 8 Pot 7 a E rived ] SCOtr == emene (24) 

where = \/(4r)— wy). 

The experimental values of the ionization potentials, taken from the 
compilation of Bacher and Goudsmit (1932), together with the calculated 
values of 7), are shown in table 1. The values of the ionic radius 7;, as 
calculated by Pauling (1927), are also given, and it is seen that 7) is less 
than 7; only in the case of mercury; the present method is therefore 
inapplicable to this metal, and it will be omitted from further consideration. 


T9 


Table 1 


Be Mg Ca Sr Ba Zn Cd Hg 
I, 1-332 1-100 0-869 0-807 0-732 1-315 1-238 1-373 
Tan 196) 256° 3:49 3:83 4:33 1:99 217 1-87 
Fag oo ee ose 1-876 2-14 19755) 140 1-83. - 2-08 


I, (Rydberg units) is the second ionization potential of a free atom, from 
Bacher and Goudsmit (1932); 7) (Bohr units) is the calculated atomic radius 
‘corresponding to the minimum value of «(7,); 7; (Bohr units) is the ionic © 
radius, from Pauling (1927). 


§3. Tue Torat ENERGY oF THE METAL 


The eigenvalue «,(7,) is only one of several terms comprising the total 
energy of the metal, and the remaining terms will be evaluated tentatively 
as for free electrons. ‘This is, of course, a more doubtful approximation 


332 _ §. Raimes on a Calculation of the Cohesive Energves 


for the divalent metals than for the alkali metals, but previous work on 
magnesium indicates that it may not be unreasonable ; in any case 
comparison of the results with experiment might give some information 
about the true magnitudes of the energy terms, particularly the Fermi 
energy, in relation to the free-electron values. 

Using the same method as for magnesium (Raimes 1950, pp. 575-6) 
the total energy per electron of the metal is found to be 


ema 3s: 2:4 3°51 1-154 0-726 
fe) Oe | es pees petite TS A pee Sy) FE 
E(r,) ( =) a , + r2 rs r,+6-43’ a (25) 


s s, 3 s 
the terms having been kept separate for identification. The first term 
is ¢,(r,), and the last three terms, in order, are the Fermi, exchange, and 
correlation energies, all calculated according to the free-electron approxi- 
mation. As far as these terms are concerned the expression is the same, 
apart from the actual value of <9, as that obtained for magnesium (eqn. (17) 
in the above-mentioned paper). However, owing to the use of the ion-core 
field only in eqn. (1), which results in a different value for <9, an additional 
term is required here, and this turns out to be one half the self-potential 
energy of the valence-electron charge distribution within an atomic 
sphere. Assuming a uniform charge distribution, as for free-electrons, 
and remembering that there are two valence electrons to each atomic 
sphere, so that the electron density is c=6/477,?, this additional energy 
becomes, in our present units, 


1677202 (7s 2-4 
i Py eT 
0 "s 

which is the second term in eqn. (25). 

The atomic radius p corresponding to the minimum of H(r,), and the 
compressibility «, can be obtained directly from (25). The cohesive 
energy S is given by the difference between H(p) and the experimental 
value of the mean energy per valence electron in a free atom, the latter 
being minus one half the sum of the first two ionization potentials. 

Table 2 shows the calculated and observed values of p, x, S and H(p) 
for the divalent metals. The set of results marked (a) for beryllium 
were calculated as above, while those marked (b) were obtained using the 
Fermi energy calculated by Herring and Hill (1940) for this metal instead 
of the free-electron formula ; this will be discussed in § 5. 


§4. THe RELATION BETWEEN PRESSURE AND VOLUME 


Simce there are two valence electrons to each atomic cell the total 


energy per atomic volume is 2H(r,). The pressure p at the absolute 
zero of temperature is therefore given by 


ee oe 1 dE 
abe Rs cir Sag We oe ace Soe 2G) 


and this can be found directly from eqn. (25). 
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The compression is defined as 
Av w—v TNS 
—— = 2 = ]— (2) PRs Lee ate pein (27) 


where v is the volume corresponding to zero pressure. Bridgman (1941, 
1948 a) has found experimentally the compressions of the divalent metals 
for a range of pressures up to 100 000 kg/cm?, and his results are compared 
with those of the present theory in figs. 1 and 2, the full lines being the 
theoretical curves and the broken lines the experimental curves. 
Beryllium is omitted owing to the unreliability of both the theoretical 
and experimental results for this metal. 


Table 2 


S —E(p) 
(kcal/mole) (kcal/mole) 
Cale. Obs. 
528 556 
429 461 
401 431 


366 = 3398 
614 655 
584 622 
619 


707 
708 


p is the atomic radius corresponding to the minimum of E(r,); « is the 
compressibility ; S is the cohesive energy ; and H(p) is the total energy of the 
valence electrons in the metal. The observed values are at room temperature. 

{ 

Bridgman (1948 a) has given values for magnesium up to 100.000 kg/cm?, 
but only the first half of the resulting curve is shown in fig. 2, on a larger 
scale than the others, together with the theoretical curve. The discrepancy 
between the two curves is greater than one would expect for this metal, 
which, of the divalent metals, is probably best suited to the present theory. 
However, in a later paper Bridgman (1948 b) gives the results of more 
refined measurements on magnesium, up to 30000 kg/cm*, and these, 
shown by crosses in fig. 2, agree almost exactly with the theoretical results. 
It seems likely, therefore, that at higher pressures the experimental 
curve should follow the theoretical curve more closely than does the one 
shown. . 

It should be noted that the experimental results are at room-temperature 
while the calculated ones are at absolute zero. An accurate correction 
to absolute zero over such a large pressure range would hardly be feasible 
with the experimental data available, but a rough calculation on the lines 
of that given by Bardeen (1938) shows that the correction would decrease 
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Fig, 1 
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Curves showing compression as a function of pressure. The full curves are 
theoretical, and the broken curves are from the experimental results of 
Bridgman (1948 a). The separation of the curves into two groups has no 
significance, but is merely intended to avoid confusion. 
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the ordinates of the experimental curves by less than 3% ‘ Such a 
correction would alter very little the correspondence between the observed 
and the calculated curves ; that for the alkaline earth metals would be 
slightly worsened, and that for Zn and Cd slightly improved. 

The volume v used in obtaining the calculated curves was that 
corresponding to the calculated atomic radius p, so that these curves are 
entirely theoretical—a fact which should be borne in mind when comparing 


Fig. 2 


O:l2 


0:10 


0-04 


0:02 


(Oa) 


e) | 2 3 4 

px 10+ (kg/cm?). 

Curves showing compression as a function of pressure for magnesium. The full 
curve is theoretical, and the broken curve is from the experimental results 
of Bridgman (1948 a). The crosses are more accurate experimental points 
from a later paper of Bridgman (1948 b) ; these do not go beyond 30 000 


kg/em?. 
the present results with those of Bardeen for the alkali metals. The 
discontinuities in the experimental curves for Ba and Sr at about 
20 000 kg/cm? are almost certainly caused by polymorphic transitions to 
close-packed structures ; Bridgman (1941) states that a similar transition 
for Ca has been smoothed over in his results. The discontinuities at higher 
‘pressures are presumably due to transitions of the type discussed by 


Sternheimer (1950) for Cs. 
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§ 5. Furruer Discussion OF RESULTS © 


Table 2 shows that the calculated values of p and « are in fairly good 
agreement with the observed values. The agreement is poor in the case 
of S, but this is understandable since the calculated cohesive energy is 
the small difference between two large energies, one calculated and one 
empirical. It is for this reason that the values of H(p) have also been 
tabulated, and the correspondence between theory and experiment is here 
seen to be quite good—it is clear that this is not a necessary consequence 
of the use of the empirical ionization potential J,, since H(r,) contains 
four terms of similar size, only one of which depends, indirectly, upon J,. 

The principal concern of the work has been to calculate the eigenvalue 
€, belonging to yo, and the assumptions made in doing this were essentially 
those made by Frohlich (1937) in dealing with the alkali metals, plus the 
additional one that J, can, for polyvalent as well as for monovalent metals, 
be calculated using the ion-core field only. The work on magnesium 
(Raimes 1950) indicates that these assumptions are reasonable, and it is 
believed that the values of ¢, so obtained are not greatly in error. Owing 
to the further approximations involved in calculating #, however, the 
final results do not give a true idea of the accuracy of ¢) ; rather, assuming 
that <,) is fairly accurate, the results show the effect of using the free- 
electron approximation to evaluate the remaining terms in H£, and it is 
evident from table 2 and figs. 1 and 2 that this is by no means so bad as 
might have been expected. ‘The most serious error probably occurs in 
the Fermi energy, and here it might be noted that the calculated energies 
would all be improved if the Fermi energy were slightly less than the free- 
electron value. 

The results for Be have some bearing on the foregoing remarks. The 
use of the free-electron Fermi energy, in eqn. (25), gives very poor results 
for this metal, but when the Fermi energy calculated by Herring and Hill 
(1940, p. 150, table VIII) is used, the other terms being unchanged, a 
great improvement occurs, the cohesive energy almost exactly equalling 
the observed value. The calculated value of p is improved, and that of 
« altered in the right direction, but a little too far ; this over-correction 
of the compressibility is probably a consequence of our interpolating 
quadratically from the values of the Fermi energy for three atomic radii 
given by Herring and Hill. 

The same assumption regarding the use of the ion-core field only in 
calculating y%) would permit the application of the recent method of 
Kuhn and Van Vleck (1950) to polyvalent metals. This method, which 
the authors apply to the alkali metals, is in principle very similar to 
that of Frohlich, except that assumptions regarding yj, are obviated to 
some extent by the use of several empirical atomic term values. Apart 
from the disadvantages involved in the latter, however, the computational 
labour would be much heavier, and it is doubtful whether the significance 
of the results would be greatly enhanced. 
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It was remarked in § 1 that the present method can easily be extended 
to metals of higher valency, and it has, in fact, already been applied to 
several of these. It is found that the calculated compressibilities do not 
compare so well with experiment as in the case of the divalent metals, 
but the calculated values of p and EH(p) remain quite good. Perhaps it 
should be mentioned that the method can only be applied directly to metals 
for which , and the state of the single remaining valence electron in a 


free ion, from which all other valence electrons have been removed, are 
both s states. 

My thanks are due to Professor H. Jones for his kind advice and 
encouragement at all times. 
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SUMMARY 


In this paper the theory of equivalent orbitals is used to consider the 
electronic structure of diamond. The description in terms of equivalent 
orbitals is closely related to the usual chemical description since each 
orbital is localized and corresponds to one of the CC bonds. The Bloch 
orbitals can be expressed exactly as transforms of these equivalent 
orbitals. This leads to anew and more accurate method of determining 
the Bloch orbitals and their energies. For diamond, it is found that the 
energies can be expressed as an explicit function of the coordinates in 
reciprocal space. The results are compared with those obtained by 
Kimball, using the cell method, and the agreeu.ent is good. This method 
has the advantage of relating the valence description of diamond directly 
with the zone description. 


$1. INTRODUCTION 


THe electronic structure of diamond ha _ been considered, using the 
Wigner-Seitz-Slater method of calculating the wavefunctions, by Kimball 
(1935). He has calculated the energies of the Bloch functions in terms of 
their wave vectors k and discussed the significance of his results in deter- 
mining the properties of diamond. The object of this paper is to show how 
similar results may be obtained, with less labour and in a more general 
form, by using an equivalent orbital method proposed by Hall (1951). 
This method also enables the relation between the tetrahedrally-directed 
CC bonds and the zone structure to be clarified. 


§2. EQUIVALENT AND MOLECULAR ORBITALS 


The modern orbital theory of valency is an attempt to study the elec- 
tronic structure of molecules, in their ground states, by using an approxi- 
mate wavefunction in the form of a single determinant. The elements of 
this determinant are functions of the coordinates of one electron only, 
known as orbitals, each multiplied by a suitable spin factor. In order that 
a wavefunction of this form should be the best approximation to the true 
wavefunction, it has been shown that the orbitals must satisfy certain 


* Communicated by the Author. 
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integro-differential equations. These equations do not, however, deter- 
mine the orbitals uniquely for both the form and magnitude of a deter- 
minant are unaltered if the elements are subjected to an arbitrary unitary 
transformation among themselves. Thus we may distinguish orbitals of 
several different types which, nevertheless, lead to the same total energy, 
the same total electron density etc. Two particular kinds of orbitals, 
equivalent orbitals and molecular orbitals, have been found to have special 
significance. Equivalent orbitals are relatively well localized in a part of 
the molecule and correspond to the description of the molecule in terms of 
inner shells, lone pairs and bonds. Molecular orbitals, on the other hand, 
are spread throughout the molecule and are the appropriate orbitals to 
use in describing ionization or excitation. A fuller account of these 
orbitals and their properties has been given by Hall and Lennard-Jones 
(1951). 

This theory can be applied with very few modifications to valence 
crystals such as diamond or silicon. The molecular orbitals are the 
Bloch orbitals but the equivalent orbitals do not seem to have been 
introduced previously in the theory of solids. These equivalent orbitals 
certainly exist for valence “rystals since these crystals can be described in 
terms of the usual directed chemical valencies of the elements concerned. 
For metallic crystals this is not so and, while there may be equivalent 
orbitals of the delocalized type such as there are in conjugated molecules, 
their existence as transforms of the occupied Bloch orbitals would have to 
be proved in each instance. 


§3. THe Structure or Dramonp 


In chemical terms, diamonu' consists of a lattice of carbon atoms each 
surrounded, in tetrahedral directions, by four other carbon atoms and 
with all the interatomic distances equal. The equivalent orbital descrip- 
tion is also fairly simple since there are only two equivalent sets, one of 
inner shell orbitals around the various nuclei and the other of bond orbitals 
between each pair of neighbouring nuclei. The orbitals in these sets are 
equivalent in the sense that an inner shell orbital localized around one 
nucleus is transformed into an inner shell orbital localized around any 
other nucleus by some symmetry operation of the space group and simi- 
larly for the CC bond orbitals. The equivalent orbital description serves, 
therefore, as a translation of the chemical description into quantum 
mechanical terms. ; 

Crystallographically, the atoms of diamond he on two interpenetrating 
face-centred cubic lattices. The unit cell of this structure is a parallelo- 
piped of a quarter the volume of the cube and containing only two atoms. 
Tf the three basic lattice vectors are a,, a2, a, the two atoms in the (J, m, n) 


cell will be at the positions 
a=la,+ma,+7Nas, 
a’=(I+})a,+(m+})ast (n+4)as- 


2AZ 
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$4. THe ENERGY ZONES 


For certain purposes the descriptions of diamond already given are 
insufficient and the molecular or Bloch orbital description is required. 
This can be derived most easily from the equivalent orbital description. 
‘The method depends on the properties of the matrix e,, which measures 
the interaction between the equivalent orbitals x, and y,. The matrix 
T',, which reduces e,, to diagonal form is the matrix which transforms the 
equivalent orbitals to molecular orbitals, and the eigenvalues H,, are the 
energies of the corresponding molecular orbitals 4, This property is 
used as the definition of molecular orbitals. The problem becomes, 
therefore, the solving of the system of homogeneous linear equations 


Le pgl tg—L pT 1,=9 (r=1, Dire) s a (4.01) 


for each value of ¢. The unknowns 7',, are partly determined by the 
symmetry of the crystal and this simplifies the problem enormously. 

If the equivalent orbitals for the CC bonds of diamond are, in fact, 
localized, then the interaction e,, between two non-neighbour orbitals will 
be zero. Although these interactions may not be exactly zero, it is 
reasonable to assume that, in the neighbourhood of the equilibrium 
internuclear distance, they are of negligible importance compared with the 
remaining interactions. Consequently we shall assume that 


Cp =, AE ROE SY Cena 3) SA ee a 
e, =, if x, and x, are first neighbour bonds, . . (4.03) 
=0 otherwise. a eer ee iy Ae a 


Apart from the original assumption of a single-determinant wavefunction, 
this is the only approximation made in this method and, if required, it 
could be relaxed and higher order interactions included. There is no loss 
of generality in putting the interaction of inner shell orbitals and bond 
orbitals equal to zero so we may omit any further reference to the inner 
shell orbitals. 

In each unit cell there are four CC bonds which may be labelled 1(1, m,n), 
2(1, m,n), 3(1, m, n) and 4(1, m,n). The matrix 7',, can then be written 


Toe=Ag,. if s=1(0, m,; 

=B,,. if s=2(1,m, 
3(1, 
4(/ sm, 


n), 
n) 
s,m) 

) 


nr). 


om (4.05) 


=) 


if s= 


? 
Pe 


The translational symmetry suggests the substitutions 


(4.06) 
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where q is a vector in the reciprocal lattice, viz. 
q=/fb\-+gb.-+-Ab,, . . 1. 2... (4.07) 
and ape Die ewe ete entre ce od as (4,08) 
The effect of this substitution is to reduce the infinite set of equations 
(4.01) to the finite set 
(a—B)A,+8(1+0-#)B, + 0(1 +e), 4+ 6(1-+e-™)D,=0, 
WL +e%) Ag+ (a—B)By+B(1+e% 0, + 6(1-+e-D, =0, ‘ 
b(1+-e'”)A,+b(1+e"-¥)B,+ (a—E)C,+b(1+e%-")D=0, . ( .09) 
b(1+-e")A,+b(1+e"-#)B +5(1 +el—19)0 + (a—E)D,=0. 
The eigenvalues of the original matrix e,, are, therefore, the roots of the 
equation 
a—EH b(l+e-%) b(1+e-%)  6(1+e-%*) 
b(1+e*) a—K b(l+e%-"7) b(1+e%%) 
b(l+e'%) b(1+ets-*) a—K b(1+et9—th) 
b(1+e') b(1-+eth-tf) b(1+e%—i9) O@—H 
For each value of q, there are four roots to this equation so that the values 
of # fall into four zones. ‘Since each # is a triply-periodic function of q, 
it is only necessary to consider its form inside a single Brillouin zone. The 
first Brillouin zone for a face-centred cubic lattice has the form of a trun- 
cated octahedron. Hach root H may be considered inside this zone or the 
root of lowest energy may be considered in it and the remaining 
roots in the second, third and fourth zones outside it. 
The eqn. (4.10) has a double root at 


a= ob ee be a, Meera (4,12) 
so that the remaining roots can be found explicitly as 
H=a+b[2+,/{4+2 cos f+2 cos g+2 cos h+2 cos (f—g)+2 cos (gy —h) 

‘ | E208 (h—/) ko meters. 7. (4.12) 
This means that two of the roots are independent of q and so remain 
constant throughout their zones. The other two roots are symmetrical 
with respect to 


=0. (4.10) 


Bas 20 peur 28 es, 9 (4.13) 


and so have contours of exactly the same shape. 


§5. THE PLANE SECTION BY h=0 


The main feature of the energy surfaces inside the zones may be found by 
taking cross-sections of the zones. The section by h=0 is one of the most 
useful since, by symmetry, there are five other planes giving exactly the 
same cross-section. 

On the plane h=0 the energies are 

E=a+b[2++/{6+4 cos f+4 cos 9+ 2 cos (g—f)}], . (5.01) 
H=a—2b (twice). (5.02) 
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If we choose to measure the energy from an origin at a and in units of b, 
both of the parameters may be eliminated and the energy expressed as a 
pure number. This means that the shape of the energy surfaces is 
independent of any particular assumptions about the values of a and b. 
Eqns. (5.01) and (5.02) now become | 
&=244/{6+4 cos f+4 cosg+2cos(g—f)}, . - - (5.03) 
é=—2 (twice). <p eel ae Os 
The contours given by these equations are shown in fig. 1. The positive 


root of (5.03) is the lowest and is plotted inside. The negative root is 
plotted outside the first zone and the constant roots (5.04) outermost. 


Energy of the molecular orbitals as a function of q over the plane h—0. 


Because of the symmetry of the roots, the energy contours appear to be 
continuous across the boundaries of the first zone. The only direction in 
which they are continuous is f=g. 

A comparison of fig. 1 with fig. 3 of Kimball’s paper shows that there is 
close agreement between the two treatments as to the general shape of the 
energy contours, but that there are a number of differences in details. 
This figure is more general than Kimball’s in that it applies, not only at the 
equilibrium internuclear distance, but also on both sides of it. It ceases 
to be a good approximation only when, at large or very small internuclear 
distances, the interaction of second-neighbour equivalent orbitals is com- 
parable with the first-neighbour interactions. Again, since Silicon forms 
a valence crystal with the same structure as diamond, the equations and 
the figure are equally applicable to it. 
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§6. NumeEricaL REsuuts 


To convert the results already obtained into physical units the values of 

a and 6 at the equilibrium separation are required. Empirical values of 

these parameters have already been deduced from the ionization potentials 

of paraffin molecules (Hall 1951, Lennard-Jones and Hall 1952). They 
are, in ev 

a=—12-978,  b——1-458. eee G101) 


The energy levels range, therefore, from a+ 6b—=— 21-726 at the origin to 
a—2b=—10-062 at the outside. 


These results do not agree exactly with those of Kimball. The values 
Ga W198 8 ba 1471 ee ee (6,02) 


would give almost exact agreement with his results but he himself admits 
that his results are uncertain, to the extent of a small change of origin, 
so the a of (6.01) is to be preferred. The agreement between the two values 
of b is most satisfactory. 


§7. CONCLUSION 


It has been convenient, in this paper, to start from the equivalent 
orbital description of the crystal and to derive the Bloch or molecular 
orbital description from it by a suitable transformation. This trans- 
formation can be determined explicitly by completing the solution of the 
eqns. (4.09). It is then possible to find the inverse transformation 
and to use it to transform Bloch orbitals, obtained, for example, by the 
Wigner-Seitz method, into equivalent orbitals corresponding to the bonds. 
Thus the orbitals used by Kimball can be transformed into equivalent 
orbitals each localized around two neighbouring carbon atoms in accord- 
ance with chemical expectations. : 

The author is indebted to Sir John Lennard-Jones for suggesting this 
problem. 
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XXIX. Separation of the Boundary Layer at a Slightly Blunt Leading Edge 
in Supersonic Flow 
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SUMMARY 

Schlieren photographs have been taken to show the supersonic flow 
past a two-dimensional wedge placed at an incidence large enough to 
give an expansion region on one side. The leading edge was blunt in 
all cases, and its thickness ¢ was varied between 35 and 200 microns. 
For values of ¢ greater than about 80 microns the photographs showed 
a weak shock wave on the downstream side of the expansion region. 
A shock wave in this position indicates that there is a separation of the 
boundary layer near the leading edge, followed by re-attachment. The 
Reynolds number was about 1:3 105 per em, giving a value of about 
1000 for t=80 microns. It is concluded that separation of the boundary 
layer occurs when the Reynolds number, based on the thickness of the 
leading edge, is greater than about 1000. For smaller Reynolds numbers 
there was no evidence of separation. 


$1. INTRODUCTION 


In an earlier paper by Bardsley (1951) an account was given of the 
supersonic flow past a wedge with a nearly sharp leading edge, the wedge 
being set at an incidence large enough to give an expansion region on 
one side. It was found that the weak shock wave that occurred immediately 
upstream of the expansion region was caused by the slight bluntness of the 
leading edge, although the thickness of the leading edge in these experiments 
_was only about 8 microns. Because no shock wave could be seen on the 
downstream side of the expansion region, it was concluded that there was 
no separation of the boundary layer. 

In Bardsley’s experiments the leading edge was very thin and it was 
difficult to determine its exact shape. The enlarged shadowgraphs that 
were reproduced in the paper: showed that the edge was irregular and not 
rounded ; the flow past such an edge, at large distances from the edge, 
would probably be similar to that past the shape shown in fig. 1. 

Holder, Tomlinson and Rogers (1949) have investigated the supersonic 
flow past two-dimensional blunt-nosed plates, with leading-edge thick- 
nesses of the order of 1 em. They found that for nose sections at least as 
pa a I ee 


* Communicated by the Authors, 
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blunt as a semicircle there was a separation of the boundary layer on the 
downstream side of the ‘shoulder’. An oblique shock wave occurred 
where the boundary layer became re-attached after the separation. This 
will be called the ‘ second shock wave’ in this paper, to distinguish it 
from the bow shock-wave. 

The experiments of Holder et al. were made with the plates set at zero 
incidence, but it is likely that separation of the boundary layer would 
also have been observed if the plates had been set at such an incidence 
that an expansion occurred on one side. The experiments of Holder et al. 
were made with a leading-edge thickness about 1000 times as great 
as in Bardsley’s case. This suggests’ that the separation of the boundary 
layer on a blunt-nosed plate may depend on the Reynolds number based 
on the thickness of the leading edge. In the experiments of Holder et al. 
this Reynolds number was about 100000 and the boundary layer 
separated ; in Bardsley’s experiments the Reynolds number was about 
100 and there was no evidence of separation. (The possibility of 
separation occurring in the latter case, with a second shock wave either 
too weak or too short to be detected by the schlieren system, is considered 
in the present paper.) 

The experiments described below were made to investigate the effect of 
leading-edge thickness on the separation of the boundary layer. 


§2. NoTaTION 


d distance from surface to point at which shock wave attains 
maximum strength. 
h height of separation ‘ bubble ’ (see fig. 11). 
k constant in eqn. (4). 
1 length of separation ‘ bubble ’ in direction of stream (see fig. 11). 
M Mach number. 
- p pressure before shock wave. 
Ap rise of pressure in shock wave. 
r radius of curvature of concave part of ‘ bubble’ profile (see 
fig. 12). 
R, Reynolds number based on thickness t. 
s=Ap/p strength of shock wave. 
t thickness of leading edge (see fig. 1). 
u velocity outside boundary layer. 
x distance from leading edge. 
y distance from wedge, measured normal to undisturbed streams 
y ratio of specific heats. 
3 thickness of boundary layer (without separation). 
displacement thickness of boundary layer (without separation). 
inclination of tangent at point of inflexion to surface of wedge 
(see fig. 12). 
# micron (=0-001 mm). 
y coefficient of kinematic viscosity. 
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$3. EXPERIMENTAL DETAILS 


The experiments were made in the intermittent supersonic wind tunnel 
described by Bardsley and Mair (1951). The working section was 12-7 cm 
high and 10-2 cm wide, and the Mach number was 1-96. The stagnation 
pressure was approximitely atmospheric and the Reynolds number in 
the working section was about 1-3 10° per cm: A silica-gel drier was 
used at the intake to the tunnel, and the absolute humidity was always 
below 0-0005. 

The wedge used was 1-27 cm long in the direction of the stream, had an 
apex angle of 10-0° and extended across the whole width of the tunnel. 
It was supported by two arms on the downstream side, and the angle of 
incidence could be changed by altering the position of these arms. If the 
experiments had been made at a low incidence, with a fairly strong bow 
shock-wave on both sides of the wedge, the second shock wave occurring 
after the separation might have been obscured by the bow shock-wave. 
The experiments were therefore made at incidences high enough to give 
an expansion region above the wedge, so that in cases where the boundary 
layer separated this expansion region lay between the bow shock-wave and 
the second shock wave, and the latter could be easily seen. 


Fig. 1 


alta 


The wedge was made of high-carbon steel, ground and lapped to the 
shape shown in fig 1. For the first experiments the leading-edge thickness 
t was 0-2 mm (200). By lapping one face of the wedge, this thickness 
was then reduced in stages without changing the shape of the leading edge. 
Thus for all the values of t the shape was as shown in fig. 1. The dimension 
¢ was determined, in each case, by measuring the wedge thickness at 
several distances from the leading edge. It was found that the values 
of ¢ obtained by this method were consistent within about 15, and it is 
estimated that the possible error in ¢ is about +15 [Le 

The knife edges in the schlieren system were horizontal and were 
arranged so that increasing density in an upward direction was represented 
by a dark region. A flash tube having an effective duration of a few 
microseconds was used as a light source because it was available ; there 
was no evidence of any serious unsteadiness in the flow. 

Schlieren photographs were taken for leading-edge thicknesses of 
200, 110, 60 and 35 microns, with the upper surface of the wedge inclined 
at 5° and at 10° to the direction of the undisturbed stream. 


Section of front part of wedge. 
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To ensure that the optimum conditions for detecting the weak second 
shock wave were realized, all the photographs were taken with two 
different sensitivities of the schlieren system. It was found that the 
smaller of the two sensitivities gave better results in all cases. The 
photographs taken with this sensitivity are reproduced as figs. 2 to 9 
(Plates IX and X*), and details are given in table 1. 


Table 1 
Fig. No. 2 3 + 5 6 7 8 9 
Stream deflection De Lae lane Mg 10° Oe Te — aoe 
¢ (microns) ZANDE) SRINEOE AO satay ON Ee eh” CRE 
R, 2700 1500 800 500 2700 1500 800 500 


The flow is from left to right and the position of the back of the wedge 
is shown by the light expansion region on the lower side ; the supporting 
arms continue downstream beyond the edge of the field of view. In some 
of the photographs the rubber sealing material between the ends of the 
wedge and the glass windows obstructs the field of view slightly on the 
lower side ; care was taken to ensure that this did not occur on the upper 
side of the wedge. 


$4. INTERPRETATION OF THE PHOTOGRAPHS 


The photographs resemble those reproduced as figs. 4, 5, 7 and 8 in 
the paper by Bardsley (1951), except that a second shock wave can be 
seen on the downstream side of the expansion region for leading-edge 
thicknesses of 200 and 110 pw (figs. 2, 3, 6 and 7, Plates [IX and X). For 
thicknesses of 60 and 35 » (figs. 4, 5, 8 and 9, Plates [IX and X) no such 
second shock wave can be seen. The second shock wave extends for a 
considerable distance, even for the smaller value of t at which it is visible 
(110). Since the shock wave cannot be seen at all when ‘60 p, this 
suggests that an important change occurs in the flow as ¢ varies between 
60 and 110. It may also, be noted that the appearance of the second 
shock wave, and the value of ¢ at which it becomes visible, are the same 
for stream deflections of 5° and 10°. 

In cases where there is a separation of the boundary layer, followed 
by re-attachment, the flow is probably similar to that shown in fig. 10. 
A second shock wave is formed near the point of re-attachment, and the 
existence of this shock wave is an indication of boundary-layer separation. 
Hence it is concluded that separation occurs for values of ¢ greater than 
about 80, while for smaller values of ¢ there is no separation. 

Values of R, the Reynolds number based on the leading-edge 
thickness f, are given in table 1. The variation in flow pattern with 
the thickness ¢ must be interpreted as a change depending on this 
Reynolds number R#,. Thus the experiments show that separation of 
the boundary layer occurs for values of R, greater than about 1000. 


ee ee 
* For plates see end of issue. 
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For smaller values of R, no second shock wave can be detected, indicating 
that there is probably no separation. However, it is still necessary to 
consider the possibility that for small values of t the second shock wave 
may be too weak to be detected by the schlieren system. This question 
is discussed below. 


§5. STRENGTH OF THE SECOND SHOCK WAVE 


By consideration of the sensitivity of the schlieren system, Mair (1951) 
has shown that it should be possible to detect a shock wave having a 
strength (s=Ap/p) of order 0:01. At the Mach number of these experi- 
ments the deflection of the stream at a shock wave of this strength is 
only about 0-2°. Thus we must consider whether it is possible for separa- 
tion to occur in such a way that the deflection of the stream at the second 
shock wave is less than about 4°. 

For a preliminary consideration of the strength of the second shock 
wave, it will be assumed that there is a sharp corner at the point of 
re-attachment, as shown in fig. 11. Let the separation ‘ bubble’ have 
a length / in the direction of the stream and a height h above the surface 
of the wedge. 

Examination of figs. 4, 5, 8 and 9 (Plates [X and X), where no second 
shock wave is visible, shows that if there is a separation, / cannot be greater 
than about 2mm. (This is an upper limit, given by the apparent width 
of the expansion region at the wedge.) If there is a separation in these 
cases, the height of the bubble cannot be less than the thickness 6 of a 
normal unseparated boundary layer on a flat plate at a distance //2 from 
the leading edge. As shown by Bardsley (1951) the displacement thickness 
of a laminar boundary layer on a flat plate at a Mach number of about 
2 is given approximately by 


9, d(vajuyia? 0 ee ey 


At this Mach number the total thickness 3 is about twice the displacement 
thickness. This is shown, for example, by Cope and Hartree (1948). 
Hence 

Sim B(varlah Aen es. a ie) oie ee 


If /=2 mm we may put x=//2=1mm. Then 60-1 mm, and if A is 
not less than 6 the ratio h/l for this value of J cannot be less than 1/20. 
For any smaller value of / the minimum possible value of h/l is greater 
than this, so that the inclination of the downstream end of the bubble 
to the surface of the wedge cannot be less than 1/20 radian (i.e. about 3°) 
for any possible value of /. The strength of a shock wave giving a 
deflection of this magnitude is about 15 times that of the weakest shock 
wave that can be detected with the schlieren apparatus. 

These calculations suggest that any separation of the boundary layer 
should be accompanied by a second shock wave that can easily be detected, 
but they ignore the fact that the true shape of any separation bubble that 
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May exist is probably as shown in figs. 10 and 12, with a concave 
transition curve on the downstream side instead of a sharp corner. This 
will now be taken into account. 

For the present, the variation of entropy behind the bow shock-wave 
will be neglected. In fig. 12, the shock wave starts at the same point as 
the envelope formed by the Mach lines from the concave curve, and lies 
a little upstream of the first Mach line in the compression region, i.e. just 


Fig. 10 
FREE STREAM a 
x 
DIRECTION “y 


SEPARATION 
“BUBBLE” 


Flow past wedge with blunt leading edge. 


inside the expansion region formed by the Mach lines from the convex 
part of the bubble. Let be the angle between the upper surface of the 
wedge and the tangent drawn at the point of inflexion on the separation 
bubble. By calculating the positions of the Mach lines for a typical 
case it was shown that, if A was small, the shock wave penetrated such a 
short distance into the expansion region that the latter could only have a 
small effect on the strength of the shock wave. Thus the deflection at 
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the strongest part of the shock wave was nearly equal to A, a result that 
was also obtained by using the first-order theory given by Whitham (1951). 
This means that the previous discussion of the simplified flow shown 
in fig. 11 is applicable to the more realistic case shown in fig. 12. 


SEPARATION 
“BUBBLE” 


Simplified form of separation ‘‘ bubble”, with sharp corner at 
point of re-attachment. 


Fig. 12 


POINT OF A 
INFLEXION STOR AOR 


SHOCK WAVE 


TANGENT AT 
POINT OF INFLEXION 


e a 
Separation “ bubble ” with concave transition curve on downstream side. 


It is of interest to calculate the order of magnitude of the distance d 
from the surface of the wedge to the point at which the shock wave attains 
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its maximum strength. By calculating the Mach-line pattern for a 
typical case (still neglecting the variation of entropy behind the bow 
shock-wave) it was shown that the distance from the point at which the 
shock wave started to the point of maximum strength was small compared 
with the distance d. Hence d is approximately equal to the distance from 
the wedge to the point at which the shock wave starts. Meyer (1948) 
has shown that for a circular-are boundary of radius r this is given by 


d= | ee fi Mm~?2-2 
= eae |= Gori 2- 2) eee eee (3) 
(The presence of an expansion region on the upstream side of the 
compression region has no effect on this distance.) 

In the limiting case when the second shock wave is only just detectable, 
the concave part of the boundary has a length of order 1 mm and subtends 
an angle of order 0-01 radian, so that its radius of curvature r is of order 
100 mm and the distance d is of order 10mm. This means that the 
maximum shock-wave strength as calculated above occurs far enough 
from the wedge for the shock wave to be easily seen. 

The effect of variation of entropy behind the bow shock-wave will 
now be considered. For the present purpose we may define the * entropy 
boundary-layer ’ as the layer adjacent to the wedge in which the increase 
of entropy due to the effect of the bow shock-wave is greater than 10% of 
the maximum increase of entropy (occurring at the normal part of the 
bow shock-wave). Consideration of the photographs given by Holder 
et al. (1949) shows that the thickness of this entropy boundary-layer is 
of the order 10¢. But inthe important case where the second shock wave 
is only just detectable the distance d is at least 1007. This means that the 
thickness of the entropy boundary-layer is small compared with the 
distance from the surface to the point at which the shock wave starts. 
Thus the Mach-line curvature caused by the variation of entropy can 
only have a small effect on the strength or position of the second shock wave. 

It is interesting to note that if t=0-1 mm the thickness of the entropy 
boundary-layer is about 10 times the thickness of the ordinary ‘ viscosity 
boundary-layer ’ at a distance of 1 mm from the leading edge (assuming 
no separation). 

It may be concluded that, if the boundary layer separates, the second 
shock wave that is formed near the point of re-attachment will be strong 
enough to be detected by the schlieren system. Thus the absence of any 
visible second shock wave for R,<1000 indicates that the boundary layer 


does not then separate. 


§6. STRENGTH OF THE Bow SHOCK-WAVE 


The angle between the bow shock-wave and the incident stream was 
measured on the photographs at various distances from the leading edge. 
From these angles approximate values of the shock-wave strength were 
calculated, but because the shock-wave angles were only slightly greater 


352 Separation of the Boundary Layer at a slightly Blunt Leading Edge 


than the Mach angle the errors in these calculated strengths were rather 
large. Hence a small error in measuring the shock-wave angle produced 
a relatively large error in the calculated strength. 

The theory given by Friedrichs (1948) shows that the strength of the 
bow shock-wave at large distances from a blunt-nosed body should be 
inversely proportional to the square root of the distance y from the body, 
measured normal to the stream, so that 


Sho oy aa 5 es ee 


where k is a constant for geometrically similar bodies. Values of k were 
calculated for all the different cases, for various values of y. Although 
the results showed considerable scatter (as would be expected because of 
the large errors in determining s), there was no systematic difference in 
k between the cases with and without separation of the boundary layer. 
This result is to be expected if h (see fig. 11) is small compared with the 
distance from the leading edge to the point of maximum bubble height, 
even if h is not small compared with ¢t. The effective thickness of the 
leading edge is then only slightly influenced by the separation bubble. 


§7. CONCLUSION 


It has been shown that the boundary layer separates if the Reynolds 
number fF, is greater than about 1000, while for smaller values of R, 
there is no separation. This result disagrees with the usual conclusion 
from laminar boundary-layer theory that the position of separation 
should be unaffected by changes of Reynolds number.. The explanation 
of this is that when the surface curvature is very high (as at the shoulder) 
the assumptions of boundary-layer theory are not valid, especially at 
such low Reynolds numbers, so that the occurrence of separation can 
depend on the Reynolds number. 
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SUMMARY 


Bounds are obtained for the error committed when any approximate 
theory of plasticity is used to estimate yield-point loads. It is shown, 
for example, that, if the Mises yield function and potential are used to 
approximate the Taylor-Quinney data for copper, the error is less than 
1%, for a wide variety of loading programmes. ; 


In applications of the mathematical theory of plasticity it is usually 
necessary to assume a yield function and plastic potential which are 
convenient approximations to the available experimental data. The 
purpose of this note is to state explicit bounds to the error committed 
when estimating yield-point loads with any such approximate theory. 
We suppose that the arbitrary scale-factor in the plastic potential 
is so chosen that the corresponding surface in stress-space is inscribed 
in the yield surface. We use the maximum work principle (Hill 1948) 
for a plastic-rigid body caused to deform, in whole or part, by forces F per 
unit area prescribed over a portion S, of the surface, and velocities u 
prescribed over the remainder S,;. The principle states that 


[ F.udSp> | F*.udSy, 


where F* corresponds to any arbitrary equilibrium distribution of stress, 
‘not. violating the yield-condition anywhere in the body and such that 
F*—F on S,. This is valid if the yield function and plastic potential 
are identical at each point of the body. When they differ, the principle 
is still true provided that the arbitrary stress, represented by a vector o* 
in stress-space, is restricted to lie within the potential surface. For, if 
the actual stress-vector o cuts the potential surface at the point 
Ao(0<A<1), then 
(s—o*) . €>(\o—o*) . E50, 


since the strain-rate vector € is parallel to the outward normal at the 

point As, and the actual rate of work o .€ is positive. The result follows 

on integrating throughout the body. 

rT Le Le 
* Communicated by the Author. 
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For the given boundary conditions let W denote the rate of work of 
forces over S,, in a yield-point state associated with some given theory 
or experimental data. W is independent of the previous loading 
programme if the yield function and potential are identical (Hill 1951), ~ 
but may not be if they are different (so far as is known at present). For 
simplicity, suppose that the respective shapes of the yield and potential 
surfaces are constant throughout the body, though the degree of 
hardening may vary (as measured locally by the absolute size of the yield 
surface). For a hardness distribution similar to this given one, we can 
calculate the rate of work of forces over Sz in a yield-point state 
associated with any convenient choice of yield function and potential. 
Let W,, be the work-rate for any chosen harder material whose potential 
surface circumscribes the given yield surface, and W, the work-rate for 
any chosen softer material whose yield surface inscribes the given 
potential surface. Then, by the maximum work principle, 


Wy<W<Wy. 


In particular, the closest upper bound is obtained when the chosen 
yield function and potential are both equal to the given yield function, 
and the closest lower bound when they are both equal to the given 
potential. For certain given variations of u, or for prescribed load-ratios 
over the whole surface, the equation may be replaced by a similar one 
for the yield-point loads. It is to be noted that, when F=0 on Sz, 
W/W, is equal to the ratio of the respective hardnesses when similar 
yield and potential surfaces are used for both bounds (since there then 
exist associated, yield-point states differmg only by this factor of 
proportionality). 

As an illustration, suppose that it is required to approximate the 
Taylor-Quinney data for copper and aluminium (nominally isotropic) 
by the Mises yield function and Lévy-Mises potential. The data, plotted 
in principal-stress space, are conveniently shown in fig. 1 of the paper 
by Bishop and Hill (1951). It will be observed that the Mises circular 
cylinder, when circumscribed to the actual yield surface, touches it 
along generators corresponding to pure tension or compression; the 
actual potential surface also touches it there, while the inscribed cylinder 
touches the potential surface along generators corresponding to pure 
shear. The ratio of the radii of the two cylinders is about 1-03 for 
aluminium, and not more than 1-02 for copper. If, then, when F=0 on S,, 
we use yield stresses respectively 15% and 1% less than the actual ones, 
the error in W is less than 13% for aluminium and less than 1% for 
copper. Of course, local stresses or the deformation mode may differ 
much more, depending on the particular application. Alternatively, if 
we were to take as the standard of comparison the yield function and 
potential calculated by Bishop and Hill for a polycrystal, the ratio of 
the cylinder radii would be 1-07 and the error in W 34°. 
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Finally, suppose that the Mises yield function and potential are the 
standard of comparison, and that the former is approximated by the 
Tresca yield function. The corresponding hexagonal prisms, circum- 
scribed and inscribed to the circular cylinder, are in ratio 2/4/3~1-155 ; 
by taking the mean hexagon the error in W can therefore be limited to 
less than 8°% (when F=0 on S,). A similar conclusion holds if we use 
the hexagon with sides normal to vectors representing pure tension 
(Hill 1950). The limits of error are the same, too, if the Tresca yield 
function is also used as plastic potential; work in progress indicates 
that this facilitates the solution of axially-symmetric and other 
three-dimensional problems. ‘ 
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ABSTRACT 


Photometric methods of determining the ionization produced by a 
charged particle in passing through nuclear emulsions are discussed and a 
technique based on measuring the mean width of a track is described. 
The method has been tested by applying it to a study of the heavy 
primary component of cosmic radiation. It is concluded that the photo- 
metry of tracks produced by charged particles passing through electron — 
sensitive emulsions may be used as a sensitive measure of the specific 
ionization for values of the rate of loss of energy less than 30 Mev/g/cm?, 
and greater than 200 Mev/g/em?. Similar methods can also be applied 
for identifying multiply charged particles, if they are arrested in the 
emulsions, and if the length of their tracks exceeds 200 p. 


§1. INTRODUCTION 


One of the most important technical problems in nuclear physics is to 
measure the ionization produced by a charged particle in passing through 
matter. In nuclear emulsions and for particles of small charge, this 
ionization can usually be expressed as a function of the number of developed 
grains along the trajectory of the particle. Thus, it was shown by Fowler 
(1950) that for tracks in Ilford G5 emulsions with a grain density less than 
100 grains/100 1, the number of grains is, to a sufficient approximation, 
proportional to the energy loss of the particle. For larger values of the 
grain density, however, the effect of ‘ saturation ’ becomes important, and 
the subjective difficulties of counting make the results unreliable. This 
difficulty can sometimes be overcome, for the larger values of the specific 
ionization, by grain counting in less sensitive or under-developed emulsions. 
With emulsions of low sensitivity, however, we may expect a non-linear 
variation of grain density with energy loss so that careful calibration is 
necessary. 

In regions where grain counting becomes impossible, an alternative 
method of measuring the ionization is to count the number of 5-rays along 
the track. In electron sensitive emulsions, 5-rays of energy above 
~15 kev produce short tracks containing three or more grains, and their 
number per unit length can also be taken as a measure of the specific 
ionization of the particle, (Frier, Lofgren, Ney, Oppenheimer, Bradt and 


* Communicated by Professor C. F. Powell, F.R.S. 
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Peters 1948). This method has been of great importance in extending the 
range of values of the specific ionization which can be measured, but it 
suffers from the fact that it is sometimes difficult to decide whether a 
particular electron track is associated with that of the multiply charged 
particle. It is also difficult to maintain a constant criterion to ensure 
accepting only those $-rays with an energy above a certain lower limit. 
Finally, if the number of d-rays/100u exceeds 20, counting becomes very 
difficult and the results are subject to large errors. | 

An alternative to these methods is to make photometric measurements 
of the obscuration produced by a track, and very encouraging results have 
been obtained with electron sensitive emulsions 200 p thick by von Freisen 
and Kristianson (1950). These authors showed by such methods that it is 
possible to distinguish between 7 and » mesons stopping in the emulsion. 
Other interesting results have been obtained by Blau, Rudin and 
Lindenbaum (1951) in studies of tracks of the heavy component of the 
primary cosmic radiation, using Kodak NT4 emulsions 80 pu thick. 

The present experiments were undertaken to investigate the possibility 
of photometric methods in electron sensitive emulsions 400 p‘thick. The 
extension of the methods to such thick emulsions is important because they 
provide favourable conditions for energy measurements by the scattering 
method. 

§2. EXPERIMENTAL 
Direct Measurement of Obscuration Produced by a Track 

Measurements have been made on Ilford G5 Nuclear Research emulsions 
400 yw thick using the procedure of von Freisen and Kristianson. An image 
of the track was formed in the plane of a narrow slit of which the long 
axis was parallel to the track, by an objective with magnification x 95. 
The light passing through the slit was directed on to the cathode of a 
photo-multiplier tube and the resulting anode-current was measured by a 
galvanometer. A galvanometer reading was first taken when the image of 
the track was contained within the area of the slit. After moving the 
plate approximately 50 u in a direction perpendicular to that of the track, 
a second reading was taken to determine the background. If these two 
galvanometer readings are R, and R,, respectively, then the relation 
R=(R,—R,)/R, is a measure of the ‘ total obscuration ’ produced by the 
track. Provided that the image of the track in the slit is black, and that 
the number of background grains is small, R expresses the ratio between 
the projected area of the track contained in the slit and the area of the slit. 
This interpretation of R is valid if thin emulsions of high transparency 
are used, but owing to scattering of light, difficulties are met in experiments 
with thick emulsions. 


Effect of Scattering of Light in Thick Emulsions 

In order to examine the effect of the scattering of light, measurements of 
R were made, at different depths in the emulsions, on the tracks of particles 
which traversed it.. Fig. 1, curve F, shows the values of R for a typical 
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track, produced by an oxygen nucleus moving at a relativistic velocity, at 
different depths in the emulsion. It will be seen that the results suggest a 
marked decrease in the opacity of the track with increasing depth. It 
follows that the value of R can be taken as a measure of the projected area 
of the track only in a region very near the surface of the emulsion. Even in 


Fig. 1 


Total obscuration (A) or half-width (B) in arbitrary units. 
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The ‘total obscuration’ (curve A), and the ‘half-width’ (curve B), as a 
function of depth below the surface of the processed emulsion, for the track 
of an oxygen nucleus which traversed an Ilford G5 emulsion, 400 » thick. 


this region, however, measurements are not reliable because of processing 
irregularities. We have found that the function relating the value of 
R/Ryox With the depth varies in passing across a single plate so that a 
correction for the apparent decrease of opacity with depth is generally 
impracticable. 
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Several attempts to overcome these difficulties were made by decreasing 
the effect due to the scattering of light by the emulsion above the track. 
These included the use of red light and of polarized light. Changes were 
also made in the optical equipment such ag a decrease of the ‘ numerical 
aperture ’ of the objective, illumination by a narrow slit, and dark field 


illumination. They were unsuccessful and an alternative approach was 
therefore adopted. 


Distribution of Blackening of a Track 

A slit, having an apparent width in the object plane of approximately 
0-25 , was set parallel to the image of a track formed in the plane of the 
slit. By means of the horizontal stage motion, the plate was then moved 
in a direction perpendicular to that of the track through a distance of 
about 100. In effect, the image of the track passed across the slit, 


Distance from axis of track. 


* Blackening distributions’ for the tracks of three particles A, B and C, the 
ionization loss of which were 350, 820 and 1400 mev/g/cm?, respectively. 
The curves have been normalized to correspond to the same maximum 
amplitude. 


remaining always parallel with it. At regular intervals during this move- 
ment, galvanometer readings were made, and a plot of these as a function 
of the distance moved gives a measure of the ‘ blackening ’ of the track as a 
function of the distance from its axis. Examples of such observations 
for tracks of different specific ionization are shown in fig. 2. 

Experiments showed that the shape of a curve obtained by the above 
method is independent of the width of the slit if its image in the object 
plane is less than the mean diameter of a silver grain. The general form 
of the curves is also independent of the ionization of the particle producing 
the track. 
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Further, fig. 1, curve B, shows the result of measurements of the 
widths of the curves at half amplitude (‘ half-width ’), for the same track 
as that considered in curve F, and proves the that value of the half- 
width is independent of the depth of the track in the emulsion. The half- 
width may therefore be taken as a measure of the ionization of the particle. 


Fig. 3 


() 
Diagram of the apparatus. 


In order to decrease the time required for the measurements and to 
improve their accuracy and statistical weight, a semi-automatic apparatus 
was built which traces the blackening distribution on the screen of an 
oscilloscope. 


§ 3, APPARATUS 
Fig. 3 is a schematic diagram of the apparatus used in the latter part of 


this investigation. The basic component is a Cooke, Troughton and Simms 
M4000 microscope, equipped with a nuclear research stage, dA. An ‘event’ 
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in nuclear plate B, is illuminated by a substage condenser and its image is 
formed on a glass mirror surface, H. This image is viewed with an 
auxiliary low-power microscope, J, by means of a partially silvered glass 
surface, D. Through the auxiliary microscope one can observe a slit 
20 py wide by 4 mm long, contained in the mirror surface. This slit permits 
light from the objective to fall on a photo-multiplier tube, /, and the 
voltage output of the photo-multiplier is proportional to the intensity of 
the light incident upon it. This output voltage is amplified by linear 
d.c. amplifier, 17, which activates the vertical deflection coils, V, of the 
oscilloscope O. 

The horizontal deflection of the electron beam is controlled by the stage 
movement ; when the stage is moved by the rotation of a micrometer 
screw, the voltage output of the potentiometer, L, varies linearly with this 
motion, as is shown schematically in fig. 3. This varying potential is 
amplified by another d:c. amplifier, N, which activates the horizontal 
deflection coils H of the oscilloscope. 

In making measurements, the photographic plate is placed on the 
microscope stage in such a position that the track to be studied is aligned 
perpendicular to the direction of the motion of the stage indicated in fig. 3. 
The slit is then adjusted until it is parallel to the image of the track by 
observation through the auxiliary microscope. The oscilloscope pattern 
is now traced by moving the image of the track across the field of view. 

The resultant trace represents the output of the photo-multiplier tube as 
a function of the stage movement and has the same form as the curves of 
Ge ae ae 

§4. RESULTS AND DISCUSSION 

Using the apparatus described in the preceding section, a preliminary 
experiment was made to determine the change in the shape and dimensions 
of the patterns traced by the oscilloscope when the image of the track in 
the plane of the slit was defocused or was misaligned with the slit. It was 
found that the errors introduced into the measurements by small subjective 
inaccuracies, either of focussing or alignment, were negligible. It was found 
also that, to a first approximation, the area of the patterns for the same 
section of the track was independent of the focussing position. This in- 
dicates that the obscuration produced by a silver grain in the emulsion 
does not depend upon the focussing and that only the apparent area of the 
grain is affected. 

In order to determine the blackening distribution as a function of the 
distance from the centre of the track, measurements were made on a selec- 
tion of tracks of heavy nuclei of which the charge had been determined by 
measurements of 58-ray density and multiple scattering. (Dainton, 
Fowler and Kent 1951). At regular intervals along the track the oscillo- 
scope patterns were recorded on tracing paper and the ‘ average ’ pattern 
determined. Several of these curves, normalized to correspond both to a 
given value of the half-width and to a given maximum amplitude are shown. 
in figs. 4 and 5. Fig. 4 shows such normalized distributions for four 
tracks, A, B, C and D, of which the parent particles had approximately the 
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same velocity but values of the specific energy loss equal to 160, 680, 1100 
and 3000 mev/g/em?, respectively. Similar results are shown in fig. 5 for 
particles with specific ionization above 100 Mev/g/cm* but with very 
different velocities. 

The curves shown in fig. 4 indicate that, within the experimental 
errors, the form of the curve showing the distribution of blackening aCTOss 
the track is independent of the ionization produced by the particle, provided 
that it has an energy loss greater than 200 mev/g/cm?. This confirms the 
view that the magnitude of the half-width gives a satisfactory measure of 
the ionization of the particle. It also establishes in this range of specific 
energy loss the fact that the half-width of the distribution and the total 
obscuration are proportional to one another. 


Fig. 4 
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Distance from centre of track in arbitrary units. 


‘ Blackening distribution ’ for the tracks of four particles, A, B, C and D, each 
with approximately the same velocity but with specific energy-loss 160, 
680, 1100 and 3000 mev/g/cm?, respectively. The corresponding symbols 
are @, O, x and +., respectively. ‘The distributions have been normalized 


to correspond both to a given maximum amplitude and to a given value of 
the half-width, 


Since the obscuration produced by a grain does not depend appreciably 
on the sharpness with which it is in focus, any of the curves in figs. 4 and 5 
gives an approximate representation of the distribution of ionization in the 
neighbourhood of the track of the particle, except for the ‘saturated ’ 
region of the core. It is not practicable to deduce the distribution in range 
of the 6-rays from such curves, for this involves difficult and uncertain 
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corrections. Our results show, however, that this distribution does not 
change appreciably after the velocity exceeds 10!° cm/sec of the particle, 
even if it is in the extreme relativistic region. 

Another investigation was made in order to compare the obscuration 
caused by a track with the specific ionization of the particle producing it. 
Measurements of total obscuration and half-width were made on the tracks 
in Ilford G5 emulsion of 98 particles with a specific ionization varying from 
that of relativistic «-particles to that of iron nuclei in the region of 
maximum ionization. These particles were identified by measurements of 
multiple scattering and 6-ray density. Alternatively, identity was 
established in the case of particles which were arrested in the emulsion, 


either by range and 6-ray measurements or by secondary effects occurring 
at the end of the track. 


Fig. 5 


Blackening in arbitrary units. 
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Distance from centre of track in arbitrary units. 


The variation of the ‘ blackening distribution ’ of the track with velocity of the 
parent particle. Curves A, B, C, D and # were obtained from measure- 
ments on the track of a magnesium nucleus, at average residual ranges of 
150, 300, 500, 900 and 1500, respectively.’ The curves have been 
normalized to correspond to the same value of the half-width and to a given 
maximum amplitude. Curve Z coincides with the average curve of fig. 4, 
and corresponds to a velocity of the parent nucleus of approximately 


101° cm/sec. 


Values of the ‘total obscuration ’ and ‘ half-width ’ of these tracks are 
represented as a function of the specific energy loss of the particle in fig. 6, 
curves A, and.A, respectively. In this figure we have also plotted the 

following quantities as a function of the specific ionization loss : B, the 
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number of 8-ray tracks with more than 3 grains per 100 (Dainton, Fowler 
and Kent 1951); C, the number of grains/50 in under developed Kodak 
NTA emulsion (Bradt and Peters 1950) ; and D, the number of developed 
grains/50 4 in normally developed emulsions (Fowler 1950). 

The lower part of Curve A,, fig. 6, shows a linear relationship between 
the obscuration and energy loss. Curve A, is a plot of the half-width in 
arbitrary units as a function of the specific energy loss and shows the 
approximately linear relation between the two quantities, the saturation 
being significant for a specific energy loss less than 100 Mev/g/em?. The 
saturation effect is probably due to the fact that a very large percentage of 
the energy is expended in the projection of electrons of which the range is 
less than the mean diameter of a silver grain. For particles having a 
specific energy loss greater than 300 Mev/g/cem? however, the number of 
developed grains caused by 5-rays is larger than the number of grains in the 


Fig. 6 
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Specific energy loss in Mev/g/em?. 
Comparison of methods used to determine the specific energy loss of particles in 
photographic emulsions :— 

A,, ‘total obscuration’ in arbitrary units, for Ilford G5 emulsion. 

A,, * half-width ’ in arbitrary units for Ilford G5 emulsions. 
B, number of 6-ray tracks with more than 3 grains per 100, in 
Ilford G5 emulsion (Dainton, Fowler and Kent 1951) 
number of grains per 50 » in underdeveloped Kodak NTA emulsion 
(Bradt and Peters 1950). 
number of grains per 504 in normally developed Ilford G5 
emulsions (Fowler 1950) 
plotted as a function of the specific energy loss in Mev/g/em?2. 
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saturated core of the track. This fact, together with a small correction due 
to the superposition of grains belonging to different planes of the emulsion, 
allows us to explain the proportionality between half-width and specific 
energy loss. 

For particles with an energy loss greater than 300 Mev/g/em2? and with 
a velocity greater than 101°em/sec, the linear relation between photometric 
opacity and energy loss suggests the use of photometry as the most suitable 
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method of determining the ionization. 


For rates of energy loss between 


30 and 100 Mev/g/cm2 photometry is unsuitable because the method is too 


insensitive. 
Fig. 7 
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‘Integral half-width’ for tracks of nuclei ending in the emulsion. Curves 
A and B represent the integrated half-width, measured over ranges from the 
end of the track of 200 and 500 » respectively, and plotted as a function 


of the charge of the parent particle. 


Another investigation was made to determine the possibility of identify- 
ing, by similar methods, particles arrested in the emulsion. If the tracks 
of such particles are too short, other means of identification cannot be 
applied. Measurements were made of the half-width of the track at a 
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succession of points 50 apart, starting from the end of the track. The 
sum of the values thus obtained, up to any particular range from the end 
of the track, was called the ‘ integral half-width ’ and curves A and B of 
fig. 7 show values of this quantity for nuclei of different charge. For 
curve A values of the integral half-width have been determined by measure- 
ments on the last 200 p of the tracks of particles of different charge. Curve 
B represents the result of similar observations on 500 w of each track. 
These curves, together with the results of von Freisen and Kristianson, 
show that photometry provides a satisfactory method of identifying 
charged particles which stop in the emulsion. 


§5. CONCLUSION 


The results of the present work on the application of photographic 
methods in studying tracks of charged particles in photographic emulsions 
can be summarized as follows : 

A. It is confirmed that photometry is a very useful method for studying 
particles with a specific ionization so great that grain counting cannot be 
employed, and it can be applied in cases where the energy loss is less 
than 30 Mev/g/em?. In this region the only technique that can be used is 
the measurement of the total obscuration and, for such measurements, 
thick emulsions of low transparency are generally unsuitable. 

B. For studying particles with a specific energy loss greater than 
300 Mev/g/cm? the use of underdeveloped emulsions involves difficult 
calibration, and 6-ray counting is unreliable. For such particles, photo- 
metric methods, and in particular the measurement of the half-width can be 
successfully applied even with emulsions of low transparency. 

C. Multiply charged particles which produce short tracks terminating in 
the emulsion (Range >200,) can be identified by determining the integral 
half width, for this quantity varies rapidly with the charge of the particle. 


ACKNOWLEDGMENTS 


This work was originally undertaken in collaboration with Mr. J. Webb 
to whom we are particularly indebted. 

We have pleasure in thanking Prof. C. F. Powell, F.R.S., for his continued 
interest and encouragement ; Messrs. A. D. Dainton, D. W. Kent and 
P. H. Fowler, for supplying the identified tracks for this investigation, and 
for helpful discussion and advice; and Mr. M. A. Roberts and Mr. G. 
Coombs for considerable technical assistance. 


REFERENCES 


Brau, Rupin, and Linpenspavum, 1950, Rev. Sci. Instrum., 21, 978. 

Brapt and Peters, 1950, Phys. Rev., 80, 943. 

Darinton, Fowxer, and Kent, 1951, Phil. Mag., 42, 317. 

Fow er, 1950, Phil. Mag., 44, 169. 

sony Et ie Ney, OPPENHEIMER, Brapt, and Permrs, 1948, Phys. Rev., 
rales 

von FREISEN and Kristranson, 1950, Nature, Lond., 166, 686. 


[ 367 ] 


XXXII. The Low Energy Beta-Ray Spectrum oF elm 


By R. RicuMonp and H. Ross 
Cavendish Laboratory, Cambridge* 


[Received J anuary 3, 1952] 


ABSTRACT 


A thin magnetic lens beta-ray spectrometer in conjunction with a 
coincidence circuit employing a pair of photo-multipliers has been used 
to examine the low energy beta-ray spectrum of 17°Tm, which was 
expected to exhibit a forbidden shape. The spectrum was, however, 
found to have an allowed shape. The results may best be interpreted 
on the assumption that the transition is first forbidden with axial vector 
or tensor interaction, the latter interaction being slightly favoured. 


§1. INTRODUCTION 


THE low energy beta-ray spectrum of 1’°Tm (see fig. 1) might be expected 
to have a forbidden shape, because of its high log, ft value of 9-71 
(Feingold 1951) and because of the beta-gamma angular correlation 
reported by Novey (1950). An investigation to determine the degree 
of forbiddenness and to distinguish the type of interaction is described 
in this paper. 

Fig. 1 
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§2. EXPERIMENTAL DETAILS 


wlhe spectrum was examined by means of a thin magnetic lens beta-ray 
spectrometer and coincidences were recorded between the low energy 
beta group and the 84 kev gamma-ray. This method gives a more 
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accurate spectrum than the conventional one of subtracting the 80-90% 
high energy component from the Kurie plot of the complex spectrum. 
A photo-multiplier (E.M.I. type 6262) and an anthracene crystal were 
used to count the beta particles and by examination of spectra with 
known shapes the combination was shown to have an efficiency which 
was independent of electron energies in the region 0-15 Mev to 2-5 Mev. 
The gamma-rays were counted by a similar photo-multiplier used with 
a liquid phosphor (terphenyl dissolved in toluene) which was contained 
in a glass flask placed outside the spectrometer vacuum chamber. 
A polythene window in the wall of the spectrometer transmitted the 
gamma-rays from the source but completely absorbed the beta-particles. 
Experiments using brass absorbers showed that 60°% of the gamma-ray 
counting rate was due to the 84 kev gamma-radiation and the remaining 
40% to the K x-rays of 1°Yb. The flask containing the liquid phosphor 
subtended a solid angle of about 34 steradians at the source. This angle 
was sufficiently large to ensure that no distortion of the spectrum was 
caused by an energy-dependent beta-gamma angular correlation. No 
angular correlation is expected between the beta-particles and the K 
x-radiation (Novey 1950). 

Both multipliers were magnetically shielded and the resolving time 
of the coincidence circuit used with them was 10-%sec. In order to 
obtain maximum transmission the spectrometer was used with a 
momentum resolution of 12°, and corrections for this resolution were 
made near the end point of the spectrum. Since for an adequate 
coincidence counting rate it was necessary to use a thick source (35 mg cm~*) 
several isotopes giving known spectral shapes were examined using 
sources of this thickness in order to determine the effects of scattering 
and absorption in the source and of possible back-scattering from the 
polythene window. In no case was distortion of the spectrum observed 
for electron energies greater than 0-4 Mev. The coincidence spectra of 
*4Na and 1%Au gave allowed shapes with end points of 1-:36-+0-02 Mev 
and 0-973 -+-0-012 Mev respectively. As an additional check the spectra 
of °4Na and °°Y were examined with the beta counter alone. An allowed 
shape was again obtained for 24Na with an end point of 1-37--0-02 mev. 
The Kurie plot of the °°Y spectrum was transformed by the a (=p?+q?) 
correction factor (Konopinski 1943) into a straight line with an end point 
of 2-27+-0-02 Mey. All these results are in close agreement with accepted 
values. 

§3. EXPERIMENTAL RESULTS 


An uncorrected Kurie plot of the low energy spectrum of 17°Tm 
obtained in our experiments for energies greater than 0-4 Mev is shown 
in fig. 2, together with Kurie plots to which the first forbidden correction 
factors a, A,’ and A_’ (using the notation of Konopinski 1943) have been 
applied. The best straight line is given by the uncorrected plot or the 
A,’ correction, while the A_’ correction gives a rather worse line. The 
a correction and also a general second forbidden correction factor 
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(WoW)? are definitely ruled: out. (The A ,” correction differs only 
slightly from the A, and A ,’’ corrections and from the allowed shape 
owing to the predominance of the energy-independent term («Z/2R)?.) 
The high energy spectrum of !7°Tm was found to have an end point 
of 0-970 +0-010 Mev, in agreement with the results of Fraser (1949). 
Tne end point of the low energy spectrum should be the difference 
between the high energy end point and the energy of the gamma-ray, 
Le. 0-886+0-010 Mev. Table 1 gives the end points obtained by the 
method of least squares for the curves of fig. 2- 
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V/(N/p*F)= 10.) 


Table 1 
Correction Uncorrected Ate An a=p+¢_ 
End point . 
(Mev) 0-874-40-010 0-876-40-013  0-870-£0-015 — 0-800-+0-040 


Again the @ correction may be excluded, while the A_’ correction 
gives a result slightly further from the expected value than the 4,’ 
correction or the uncorrected shape. 
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$4. Discussion 

170Yb ig an even-even nucleus and its ground state should therefore 
have spin zero and even parity. Measurements of conversion electron 
intensities (Grant and Richmond 1949) and of the lifetime of the 84 kev 
excited level (Bell and Graham 1950) have shown that the 84 kev transition 
is electric quadrupole and consequently the excited level has spin two and 
even parity. A spin of three or more for 17°Tm would make the low 
energy beta group more favoured than the high energy group and this 
objection also applies to spin two and even parity. 

Spin two and odd parity would make both beta transitions first 
forbidden but would require the high energy transition to have an a shape. 
This seems to be excluded by the observations of Fraser (1949). By 
assuming that the high energy transition was of the allowed shape he was 
able to distinguish in his complex Kurie plot of 17°T'm a lower energy 
beta group of end point 0-900+-0-020 Mev which was a straight line. 
As this result agrees well with our observation of the low energy spectrum 
one may reverse Fraser’s argument and deduce that the high energy 
transition is in fact of the allowed shape. 

Of the remaining possibilities spin zero -may be excluded since rH 
‘odd parity this would lead to an a shape for the low energy spectrum or 
with even parity to a second forbidden shape. A spin of one and even 
parity would make both beta transitions allowed in contradiction to 
their ft values and (in the case of the low energy beta group) to the 
observed beta-gamma angular correlation. 1!’°T'm therefore has spin one 
and odd parity, as indicated in fig. 1, making both beta transitions first 
forbidden. 

The shell theory of nuclear structure does not seem to be very well 
established at present for odd—odd nuclei. However the stable isotopes 
169Tm (Z=69) and 171Yb (N=101) are both known to have spin } (Mack 
1950) and therefore 17°Tm may well have its odd proton in an S, and 
its odd neutron in a P,, state. The shell theory definitely predicts odd 
parity ; and the spin would be zero or one for this assignment of orbits. 
With the * parallel spins * assumption (Suess, Haxel and Jensen 1949, 
Nordheim 1950) the total spin would be zero, whereas the observation 
seems to indicate spin one. 

The angular correlation function W(@) is of the form 1+-R/Q cos* 6 
for all first forbidden beta transitions followed by a gamma-ray (Falkoff 
and Uhlenbeck 1950), and Novey (1950) has obtained for !7°Tm the 
correlation function 1—(0-259+-0-053) cos? 6. The beta transition 
I(—)+2(+-) may correspond to any of the nuclear matrix elements 
Ir, Jz, Joxr and B, (Konopinski 1943), but of these only fo xr would 
give the negative coefficient R/Q obtained by Novey. The matrix element 
Joxr is associated with the correction factors A,’ (tensor interaction) 
and A_' (axial vector interaction). Our results are therefore slightly 
in favour of the tensor interaction. It is interesting to note that Nov ey’s 

value of R/Q is in good agreement with the theory (Falkoff and Uhlenbeck 
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1950) assuming that the beta transition is first forbidden and obeys the 
tensor form of interaction, but it must be borne in mind that the 
theoretical value obtained from the figures and tables of Falkoff and 
Uhlenbeck (1950) is strictly correct only for the case Z=0. 


We wish to express our thanks to Mr. J. M. C. Scott for many valuable 
discussions. 
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XXXII. CORRESPONDENCE 


Anomalous Paramagnetism and Exchange Interaction in Copper 
Acetate 


By B. Bieaney and K. D. BowErRs 
Clarendon Laboratory, Oxford* 


[Received January 16, 1952] 


THE majority of copper salts of reasonable magnetic dilution obey Curie’s 
law with high precision down to very low temperatures. Guha (1951) 
found, however, that copper acetate, Cu(CH,COO),, H,O had a maximum 
in its susceptibility (y) at about room temperature, and at low tempera- 
tures y fell so rapidly that it would reach zero at about 50° k. In view of 
this paramagnetic resonance measurements at 0-3—0-8 cm! were under- 
taken on a single crystal, which showed that the spectrum is quite different 
from the normal. Soon after we had begun, Lancaster and Gordy (1951) 
reported resonance measurements on a powder showing similar behaviour, 
and suggested that a line observed at twice the normal g-value could be 
explained through interaction between neighbouring copper ions (cf. 
cobalt (Bleaney and Ingram 1951)). The results of our investigation, - 
and their interpretation, are given below. The measurements have been 
made principally at 90° k. , 
Copper acetate grows monoclinic prismatic crystals; the detailed 
crystallography has not been investigated by x-ray methods. The 
paramagnetic resonance spectrum is found to resemble that of a nickel 
Tutton salt (cf. Griffiths and Owen 1951), with two paramagnetic units of 


apparent spin S=1 in unit cell. Al] our measurements can be fitted to the ~ 
Hamiltonian 


H =DS?+H(S?—S,*)+B[9,H S82+92H 2tGyH yy] Je yy (1) 
using the values D=0-34 cm~1, H=0-01 em}, g,= 2-42, g,=9,=9,=2:08. 
Thus the z-axis (making an angle «= -+33° with the ac-plane) is an axis 
or nearly tetragonal symmetry, and the g-values are similar to those found 
in other copper salts. There is a zero field splitting of approximately 
0-34 cm~1, and the lines show a curious structure (see Photograph) 
parallel to the tetragonal axis. This collapses into a single line with no 
resolved structure perpendicular to this axis. 

These results can be satisfactorily explained if strong isotropic exchange 
forces are postulated between isolated pairs of ions whose individual 
magnetic axes are parallel to one another. The exchange force is small 
compared with the crystalline electric field acting on each ion, so that we 
may take these as orbital singlet states, each with effective spin }. The 


a! 
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ee couples these spins into a lower antiparallel (and diamagnetic) 
$s a with total spin S—0, and a higher triplet state with total spins — 1. 
which is paramagnetic. Only the spectrum of this latter state can be ob- 
served, of course, and it disappears at 20° k, where substantially all the 
pairs are in the S—0 state. 
If the exchange interaction is taken to be of the simple form —J S, . S,, 
the Hamiltonian may be written : 


wo X7-EX ALL, . 8,41, . 8:]—JS8, . 8. 8H, +28,)+(L.4+ 28.)]. 
(2) 


Here subscripts 1, 2 refer to the two ions; 2 is the spin orbit coupling 
coefficient. The terms are written in decreasing order of magnitude. 
Whether the crystalline field (X) is taken either as (cubic+tetragonal) as in 
the Tutton salts, or as purely tetragonal (corresponding to a likely four 


The hyperfine structure of copper acetate. 


planar bonding of copper to the four oxygens of the two acetate groups), 
the same orbital level will be lowest. The combined effect of the spin- 
orbit coupling and the exchange gives a singlet state at jJ, plus a triplet 
at —iJ. The triplet has a small: splitting, which, in terms of the 
Hamiltonian (1) may be written ) 


while the formulae for the g-values are the same as for the copper Tutton 
salts (Abragam and Pryce 1951). . 
Calculations of the susceptibility shows that our measurements of the 
g-values and of the angle « give close agreement with the anisotropy 
observed by Guha. From the temperature variation in his measurements 
we may obtain the value of J ; a very good fit is obtained if J is assumed to 
yary slightly with temperature, from — 230 cm-! at 80° K to —270 cm! at 
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300° x. Such a variation is not improbable, as J varies rapidly with inter- 
ionic distance, which will change slightly with temperature. This value of 
J, inserted in (3), gives D1 em~!, which is of the same order of magnitude 
as that observed experimentally. The discrepancy can be explained if J is 
different for the various orbital levels, since the splitting arises from the 
small admixture of these into the ground state through the spin-orbit 
coupling. Since both the splitting and the g-values are related to the 
crystalline field, the theory shows that both should have the same symmetry 
axis. This would not necessarily be the case with a splitting due to 
magnetic dipole-dipole interaction between neighbouring ions, for instance. 

Our assumption leads to a natural interpretation of the line structure 
shown in the photograph. If we add to (1) a term (S,. 4./,+4+8,. A . J) 
representing the usual magnetic hyperfine structure for the two ions (4 
being a tensor), each spectral line will be split into several components 
shifted by A(m,++m,) from the centre. Here m,, ms are the nuclear 
magnetic quantum numbers, which each range from + 3/2 to —3/2 for 
copper. It is easily seen that this gives seven components with relative 
intensity 1:2:3:4:3:2:1, which corresponds closely to the observed 
structure. Parallel to the tetragonal axis, a is approximately 0-008 cm= 1! 
while perpendicular it is less than 0-001 cm~!. These values are not 
greatly different from those in other copper salts. 

The mechanism we have postulated is akin to that of antiferromagnetism, 
except that only pairs of ions are involved, rather than a co-operative 
effect between many ions. In nuclear paramagnetism it may be likened to 
the case of ortho- and para-hydrogen. It appears to give a satisfactory 
explanation of both the fine structure and the hyperfine structure of the 
paramagnetic resonance spectrum, and of the anomalous variation of the 
susceptibility. 
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The L Absorption Spectrum of Nickel. An Anomaly 


By Y. CaucuHois 


Laboratoire de Chimie Physique, Paris* 


[Received January 21, 1952] 


THs L absorption spectrum of nickel had not been measured previously. 
In this laboratory we have obtained it with the help of a bent-crystal 
vacuum spectrograph (Cauchois 1945, 1952) and nickel foils obtained by 
beating. A gypsum crystal provides, in the experimental device, a dis- 
persion of | ev/mm. By registering on the same film, one above the other, 
the absorption spectrum and the emission spectrum from the same sample 
when bombarded by electrons, it appears that the L,,; absorption edge does 
not coincide with the emission edge. The absorption edge is displaced 
towards lower energies as compared to the emission edge. Discrepancies 
are of the order of several times 0-1 ev. 

It is well-known that emission edges and absorption. edges of a metal are 
expected to coincide and in fact do whenever experimental tests could be 
made (Cauchois 1948). Ni Ly,; seems to be the first and only case when a 
definite difference is observed. 

The measurements of absorption were made on a foil at ordinary 
temperature, while for enaission the foil was heated locally by electron 
bombardment. This might result in a broadening of the emission edge, 
together with a mean displacement (due to thermal expansion) towards 
low frequencies; this is in the opposite direction to the experimental 
effect reported here. 

The observed displacement in our view might well be connected with the 
magnetic properties of nickel. Further experiments were therefore made 
with absorbing foils which were heated somewhat above the Curie point. 
The absorption edge then comes out at higher frequencies than when the 
metal is ‘ cold’; it is practically at the same position as the emission edge. 

The sample used is of somewhat impure nickel whose characteristic data 
are not precisely known. Its expansion coefficient « may well show some 
anomalous behaviour at the Curie point, such as is known for some nickel 
alloys. A sudden decrease in « above the Curie point would displace the 
edge towards the high energies; but one would expect the order of 
magnitude of this displacement to be much smaller than the observed one. 

Further experiments are in progress. Tentative explanations will be 
proposed later. 
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XXXIV. Notices of New Books and Periodicals received 


Theoretical Physics. By G. Joos. Translated by I. M. Freman. 2nd Edition, 
(Blackie and Sons, Ltd.) Price 50s. net. 

A New and up-to-date edition of Professor Joos’s well-known book on theoretical 

physics is very welcome. There have been extensive revisions in the sections on 

wave mechanics and atomic physics generally, and a new addition of interest is a 


chapter on various technical problems. It forms a valuable textbook for any — 


student of experimental or theoretical physics, and it is convenient to find 
mechanics, vibrations, optics, electricity and quantum mechanics treated in the 
same volume and with consistent notation. , N. F. M. 


Ferromagnetisme et Antiferromagnetisme. Colloques Internationaux du Centre 
National de la Recherch Scientifique. 


Tuts is the report of the conference on magnetism held at Grenoble in July 1950, 
under the direction of Professor Neel. | Among other interesting papers the 
reviewer was particularly impressed by the review of antiferromagnetism by 
Professor Van Vleck. . N. F. M. 


Modern Magnetism. By L. F. Batrss, F.R.S. [Pp. 506.] Third Edition. 

(Cambridge University Press.) Price 30s. net. 

Iv is a tribute to this excellent work which first appeared in 1939 that it should 
already have reached the stage of a third edition. The clarity of its style and its 
deliberate, though not undue, emphasis upon the experimental aspects of the 
subject have appealed to a wide circle of readers. 

Though the second edition in 1948 attempted torbring the work up to date by 
the inclusion of new chapters, it has now been possible to accomplish this far 
more satisfactorily by recast and rearrangement in appropriate places within the 
twelve chapters of the book. Many will find special interest in the new material 
forming much of chapters 10 and 12, where the phenomenon of paramagnetic 
resonance and the recent development of the domain concept are adequately 
described ; and in the final section of the volume dealing with the author’s 
own work upon the hysteresis cycle. AS Mya 


The Preparation of Programs for an Electronic Digital Computer. By M. V. 
Witkgs, D. J. WHEELER and 8. Gitt. |[Pp. 167.] (Addison-Wesley Press, 
Inc.) : . 

Wuitst this book is designed primarily for those directly interested in the 

operation or construction of electronic computing machines, Parts I. and IT. 

may be profitably read to acquire a general knowledge both of the operation 
of such machines and, as the authors suggest, the applicability of electronic 
computers to specific fields. The book, which is confined almost exclusively 
to the KDSAC, is divided into three parts: Part I. contains a detailed 
discussion of the framing of orders and their setting in the machine, the design 
of programmes, and the use of standard ‘sub-routines’ enabling particular 
projects to be built up of previously tabulated operations. Part I. gives the 
specifications of all the sub-noutines used in the Cambridge University 


Mathematical Laboratory, and in Part TI. the programme sheets for about 


one-half of the specified sub-routines are given in full. “The second part of the 
book provides a concise and comprehensive list of the types of operation which 
may readily be undertaken by electronic machines. J. F. W. B. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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